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20  abstract  fCnntlmie  on  res-erne  %ide  if  nere^hary  and  identif\  hy  hlotk  nnrrher) 

Part  1 — We  have  surveyed  aspects  of  striation  decay  due  to  both  ionization 
falling  (and  the  subsequent  recombination)  and  diffusion  perpendi- 
cular to  the  magnetic  field.  Among  the  phenomena  assesseed  are: 

1.  Ionized  material  falling  due  to  motion  parallel  to  the  magnetic 
field,  B,  which  is  allowed  to  make  an  arbitrary  angle  with  the 
vert  leal . 
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4.  TITI.E  (and  Subtitle')  (Continued) 

Part  I — Velocity  Sliear,  the  E x B Instability  and  tlie  k ^ Power  Density 
Spect  rum. 

20.  ABSTRACT  (Continued) 

2.  Striation  falling  due  to  motion  perpendicular  to  B.  (This  occurs 
for  a nonvertical  magnetic  field  through  the  setting  up  of  electric 
fields  perpendicular  to  B dut>  to  components  of  the  neutral  velocity 
or  gravity  perpendicular  to  B.) 

3.  .'■'iffusion  perpendicular  to  B due  to  drift-dissipation  mode  tur- 
bulence . 

4.  Distance  scales  for  striation  structure  for  model  situations  of 

striation  pinching  and  striation  tip  steepening.  (The  scales  are  ! 
calculated  as  functions  of  initial  dimensions  and  diffusion  I 

coef  f ic ients . { 

I 

5.  Time  scales  for  striation  decay  due  to  the  above  processes.  j 

The  effect  of  ionization  falling  appears  to  be  capable  of  removing  * 

striation  inhomogeneities  on  a time  scale  of  order  1 hour;  diffusion 
perpendicular  to  the  magnetic  field  is  strongly  dependent  on  striation  ' 
dimensions  and  diffusion  mechanisms,  but  significant  turbulent  diffu- 
sion appears  possible  in  some  cases  on  time  scales  of  the  order  of  ten 
minutes.  There  is  some  indication  in  the  experimental  data  of  plasma 
turbulent  diffusion.  These  results  are  of  the  nature  of  estimates  i 

rather  than  detailed  calculations  and  hence  are  suggestive  rather  than  I 
def  init ive. 

Part  2 — In  this  paper  we  have  attempted  to  present  a basis  for  the  physical 

understanding  of  drift  modes  as  well  as  background  information  for  • 

their  application  to  striation  decay  problems.  If  one  neglects  the  | 

effect  of  a component  of  the  ambient  electric  field,  E,  in  the  direc-  I 
tion  of  the  background  density  gradient,  Vn  (with  botli  perpendicular  tc. 
to  the  ambient  magnetic  field,  B) , drift  modes  appear  capable  o«"  pro-  | 
venting  elongated  striatlons  from  getting  thinner  than  0.1  km  in  many 
cases  of  interest.  The  inclusion  of  electric  field  effects  for  drift 
modes  is  an  ongoing  problem. 

Part  3 — The  E x B Instability  Is  analyzecj  in  the  presence  of  plasma  velocity 
shear,  with  the  magnetic  field,  B,  in  the  z-dlrection,  the  background 
density  gradient  in  the  x-direction,  and  electric  field  components 
E (x)  and  E with  E /E  non-zero.  It  is  shown  that  the  structure  of 
E^x  B modes^under  tfies?  circumstances  can  account  naturally  for  a k~^ 
ionospheric  power  density  spectrum,  provided  one  assumes  that  the  modes 
grow  until  the  modal  density  gradient  in  the  direction  of  the  ambient 
density  gradient  becomes  comparable  with  the  ambient  density  gradient. 
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SUMMARY 


JAYCOR  work  on  striation  decay  during  the  past  year  is 

summarized  in  the  three  accompanying  reports,  "Aspects  of 

Striation  Decay,"  "A  Synopsis  of  Drift  Mode  Behavior,"  and 

_o 

"Velocity  Shear,  the  E x B Instability  and  the  k Power 
LK'nsity  Spectrum." 

The  first  report  presents  estimates  that  plasma  stri- 
ations  could  decay  due  to  ionization  falling  (and  the  result- 
ant recombination)  on  the  time  scale  of  an  hour  and  due  to 
plasma  turbulence  from  drift-dissipative  modes  on  time  scales 
of  the  order  of  ten  minutes  (under  appropriate  ambient  condi- 
tions) ; further  the  competition  between  ordinary  collisional 
diffusion  and  drift-dissipative  diffusion  is  assessed  and  con- 
ditions for  dominance  of  the  turbulent  effects  (sometimes  by 
a factor  of  100)  are  established;  in  addition  models  are  intro- 
duced to  study  the  effect  of  sensitivity  of  striation  decay 
time  to  variable  diffusion  coefficient.  Finally  it  is  estab- 
lished that  0. 1-0.2  km  structures  associated  with  anomalous 
gigahertz  scintillation  can  be  explained  by  a balance  between 
drift-dissipative  diffusion  and  convection  whereas  they  cannot 
be  explained  by  a balance  between  collisional  diffusion  and 
convection.  Hence  we  consider  that  the  first  report  serves 
as  a justification  for  the  relevance  of  plasma  turbulence  be- 
havior to  striation  decay.  The  use  of  drift-dissipative  turbu- 
lence is  to  some  extent  illustrative  and  tends  to  establish  a 
lower  bound  for  turbulent  effects. 

The  second  report  is  an  introduction  to  drift  wave 
phenomenology  with  relevance  to  the  striation  decay  problem. 

It  is  meant  to  indicate  how  properties  of  drift  modes  used  in 

the  first  report  were  obtained  as  well  as  to  indicate  a physical  I 

1 

j 


picture  for  the  instability  mechanism  and  a justification  for 
usage  of  the  dri f t -dissi pat i ve  mode  instead  of  other  drift 
modes.  It  also  indicates  some  of  the  caveat.s  necessary  when 
dealing  with  applications  of  drift  waves  to  the  striation  decay 
problem . 

The  third  report  presents  a theory  for  establishment  of 

—2 

the  k power  density  spectrum  on  the  basis  of  E B eigenmodes 

with  plasma  velocity  shear  (E  = electric  field,  n = ionized 

density,  E • Vn  / 0) . We  feel  that  this  study  of  yet  another 

form  of  possible  plasma  turbulence  is  significant  because  of 

its  direct  connection  to  the  experimentally  observed  spectrum 

and  the  resulting  possibility  of  experimental  examination,  as 

well  as  its  eventual  application  to  striation  decay.  To  our 

knowledge  although  there  has  been  much  previous  work  on  E x B 

turbulence,  this  is  the  first  work  entirely  within  the  scope 

_2 

of  E X B turbulence  which  provides  an  explanation  for  the  k 
power  density  spectrum. 
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1.  INTRODUCTION 


The  puri)()S("  of  this  i-cport  is  t,o  jiresf-nt  various  factors 
atfectinp-  striation  cifMUiy  as  well  as  to  indicate  their  sip;nifi- 
canc'O  in  ti'rms  of  possible  time  scales  for  strOation  decay  and 
spatial  structuring  of  striations.  In  the  second  st'ction 
mec'hanisnis  affiuMintt  sinkitift  of  ionization  both  parallel  and 
perpendicular  to  I^  are  present(‘d  and  time  scales  estimated.  In 
t h{?  third  section  diffusion  due  to  one  typ<'  of  plasma  turbulcuice 
is  presented  as  w('ll  as  ordinary  collisional  diffusion.  Some 
aspects  f)f  the  interj)lay  of  ditfusi<jn  with  spatial  strucOure 
and  striation  evolution  ar<?  shown.  These  are  use.'d  to  ass(>ss 
the  de|)endence  of  decay  on  diffusion  ccjeffic.ient  as  well  as  to 
infer  th<’  i^ossibility  of  plasma  turbulenc'e  diffusion  in  observa- 
tional situations.  In  the  fourth  section  our  results  are  sum- 
marizf'd  both  as  to  decay  and  spatial  structuring  of  striations. 
In  the  fifth  .secti(;n,  suM,f!:e.st  ions  for  ful  ur(>  work  are  made?. 
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m'THANlSMS  AFFE('TING  SINKING  OF  lONI'/ATION 


Coti t r i but  i ons  2.1  throuf^h  2.4  refer  to  sinkinp  ;ind  decay 
du<'  to  ionization  motion  paralliM  to  . Contribution  2 . .5  refers 
to  sinkitiK  ef  density  eniiancement  s due  to  forc’.e  components 
pe  rfxuidi  cu  lar  to  f1 . 


2.1  VERTICAL  CONTRIBUTIONS  DUE  TO  FORCE  COMPONENTS 

PARALLEI,  TO  ^ 

Th<'  separate  contributions  in  this  category  are  labelled 
V.J  , and  . 'Ve  have: 

a)  ~ ~ '■’n  D sin  D (1) 

whe  re 

= northward  comporu^nt  of  neutral  wind, 

I)  = dip  angle , 

i.s  contribution  from  vf’rtical  ct^mp.nent  of  tlu?  horizontal 
n«;utral  wind  proje,>cted  onto  an  (inclined)  magmatic  field. 

Fig.  la  illustrates  . 


b) 


v„  = 


- sin^D 

n . 

1 


V . m . 
in  1 


j 

3z 


n . (T 
1 


+ T.  ) 


(2) 


v'2  is  contribution  from  diffusion  parallel  to  B projected  verti- 
cally, n.  = ion  density,  = ion-neutrai  collision  frequency, 

k = Boltzmann's  constant,  m^  = ion  mass,  z = vertical  coordinate', 
T^  = electron  tempe;  nature , T.  = ion  ttimpf'rature , z is  vertical 
coordinate . 


c) 


g sin  1) 


^n 


(3) 


is  the  fall  from  gravity,  g = gravitational  acceleration. 


2.2  RECOMBINATION  (’IIANIBM  - GIAROE  EXHIANGK  FOLLOWED  BY 
D I SSO(’  I AT  I VE  RI'ICOMB  I NAT  I ON  ^ 1 ) 

The  loss  ratf'  is  ol  t ho  form 

- k.^  nfOg)  + n(N^)  n.  . (4) 

Here  k^  and  k^  are  constants,  n(0,^)  and  n(N^)  are  molecular 
oxyt'en  and  n i t rf)gen  densities. 

2.3  DISCUSSION  OF  CONTRIBUTIONS 

Gravity  pulls  ionized  particles  down  the  (ield  line,  at 
nipht  the  neutral  wind  pushes  ionizt.'d  particles  up  the  field 
lines  ( ■ 0),  during  the  day  the  neutral  wind  pulls  ionized 
particles  down  the  field  lines  ( > 0),  thermal  diffusion 

(proportional  to  the  ionized  particle  pressure  gradient)  can 
act  in  both  directions. 

The  ion  continuity  euuation  is: 

TT  ^ ^ ^ ^2  ^ “ [*'i  "^"2)  ^ ‘^2  "<^'2^]  "i 


wh(?re  Q represents  creation  and  z is  the  vertical  coordinate. 
During  days  rt.'lative  striation  densities  will  be  dropped  down 
by  Q as  well  as  all  other  factors  except  tliermal  diffusion.  At 
night  Q is  negligible,  Eci . (5)  is  linear  in  n.  (neglecting 
entrainment  of  neutrals  due  to  very  ti  i gh  ionized  densities)  and 


decay  time  .sc-ale.s  of  2 hours' 
rea.son able  . 


(e- folding  in  2 hour.s)  appears 


2.4  ADDITIONAL  CONTRIBUTIONS  UNDER  DISTURBED  NUCIJ-IAR 
CONDITIONS 

For  disturbed  situations  vc?rtical  neutral  winds  initially 
upwards  but  eventually  downwards  ( F.  Fa.jen  and  R.  W.  Kilb, 
private  communication)  could  play  a profound  role  in  pulling 
ionized  particles  down  the  field  lines.  Ionized  material  re- 
turning from  the  magnetosphere'  could  perform  the  same'  functie^n. 
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There  could  be  s(  roiig  et  fects  due  to  disturbed  neutral  wind 
inhomo^^eiu^i  t i<>s . 

2.5  CONTRIBUTIONS  DUE  TO  FORCE  COMPONENTS  PERPENDICULAR  TO 

These  can  involve  neutral  winds,  f^ravity  and  electric 
fields.  Thermal  diffusion  B should  be  insi  pn  i f i cant  compared 

to  that  II  B.  In  the  penerally  appropriate  limit, 

V V . 

« .-lil  « 1 . 

u)  w ■ 

C(>  cl 

(v  = e 1 ect  ron-neut  ra  1 collision  freciuency,  uj  = electron 
en  c(' 

cyclotron  freciuency)  electrons  can  only  move  with  velocity, 

cE  X B 


who  re 


E = E_^„  - 


(Here  E|  is  the  spatially  uniform,  ambient  electric  field  ('B) 

I 

and  -V I i;/  is  the  s(0  f-con.s  i .sten  t induced  electric  field  B.) 
i 

Hence,  by  quas  i -iieut  ra  1 i t y of  ions  and  eD'ctrons,  the  problem 
of  determining  ionization  density  changes  [ B is  equivalent  to 
finding  iIj.  We  take  ijj  to  be  given  by: 


E + 
o 


V X B 


where,  with  z denoting  location  along  a magnetic  field  line. 


y.  = f dz  n. 


V . e c 
in 

o)  . B 

Cl 


N = r dz  n . 


and  V is  a mean  value  for  the  field  line  integrated  neutral 
wind. 


%n  ^ 
"ci  » 
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This  o(]uation  is  based  on  Fai . ( Ki)  (jI  Flo  f<?rence  .'5  wi  Ui  the 
assumpi  i oiis  ; 


and 


1 . 


Neglect  of  hiK^tii’  order  terms  in 


2.  , V and  g are  constant  in  space  (n.b.  , v is  a 
field  line  average). 

3.  On  defining  ^ as  the  maximum  of  | j , |v|B/c 


and  ! glN^./(  Ew^.  ) , 


2T 


e E 


e f f 


s ^ 2T 
n e E 


off 


where  h^  is  the  scale  height  for  neutral  variation. 

( 4 ) 

Using  a waterbag  model , 


II 

< P 

o 

II 

> ^ ^o 

N = N.^  , p 

< p . 
() 

: ^ = ^2 

■ =*  ‘V, 

where  p is  the  cylindrical 
dir(?ction  of  B/|B|  , we  have 

radius  around 
( se(’  Ap))endix 

th€?  axis  in  the 
A for  details) 

cEl  X B 2Z., 

J.o  2 

"1  ^ 

b2  . E, 

2 

c ("-2  - 

( V X 

13)  X B , ^^1 

^2^  e(g  X B)  X B 

y2  Ej  + i-2 

- 

c E^  - 

(7) 

I f there  is  no 

inhomogeneity  in  E2 

= E^  , Nj  = N2  and 

''1 

c E|^  X B 
b2 

Hence  there  is  a vertical  velocity 


I -7 


•f 


1 


i 


c 


B 


e 

e 

— cos  D 


where  e is  the  unit  vector  in  the  east  direction.  Indications 

p 

in  Reference  2 are  that  this  is  not  a very  significant  effect. 


For  ^^‘'2  ^^^2’ 


the  vertical  component  of  vi  is 


Vi  =v  sin  I cos  I 
n 


2 . 
g cos  I 

(v  . ) 

' in/ 


(8) 


where 


The  neutral  wind  driven  contribution  to  Eq . (8)  is  shown  in 
Fip;.  lb.  It  is  interesting  to  note  on  combining  Eqs.  (l)-(3) 
and  (8)  that  in  this  limit  the  cloud  has  a net  upward  velocity 
(neglecting  altitude  differences  in  so  that  «=  ' 


- sin 
n . 


^n 


m 


7 ^ ["lOe  + T;)]  - 


\ 1 n / 


i.e.,  horizontal  neutral  winds  do  not  lower  the  cloud  (illus- 
trated in  Fig.  1)  and  gravity  lowers  it  as  if  the  magnetic  field 

were  nonexistent.  This  would  characteristically  result  in  decay 

2 

time  states  whicn  are  at  least  a factor  of  sin  D lower  than  the 
spatially  homogeneous  decay  time  scale. 

This  limit  is  of  course  suggestive  rather  thiui  fully 
realistic.  However  through  Eq.  (8),  one  sees  that  at  night 
(since  v < 0)  the  effect  of  finite  dimensions  perpendicular 
to  B is  to  increase  the  sinking  rate  of  density  enhancements 
relative  to  the  ambient. 
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*4 


Vj2,  due  to  fnrc('  components  |1  B (if  striation  conductivity 
larfte  compared  to  ambient);  horizontal  neutral  wind. 


Note  v.^  + v-2  = 0. 


Fi sure  1 . 


Illustrative  Vertical  Ionized 
Velocity  Ccjmponent.s. 


3. 


CONTRIBUTION  TO  STRIATION  DECAY  FROM 
DIFFUSION  PERPENDICULAR  TO 


The  contribution  to  strintion  decay  from  diffusion  per- 
pendicular to  l5  depends  on  two  factors; 

1.  What  are  the  mechanisms  responsible  for  diffusion. 

2.  What  is  the  striation  geometry  (shapes  and  dimen- 
sions) and  what  is  the  non - di f f us i ve  flow  on  which 
diffusions  is  superimposed. 


3.1  DIFFUSION  MECHANISMS 


Classical  ion -electron  collisional  diffusion  for  circular 
striations  perpendicular  to  l5  has  been  investigated  in  detail 
by  Kilb  and  Stagat.^^^  For  the  electron  density  n^  = n.  = 10^  cm 
and  a striation  width  of  1 km,  the  time  for  a 2.25-fold  decay  of 
density  (due  to  a cylindrical  radius  increase  of  a factor  of 
1.5)  is  about  1 hour.  For  lower  densities  and  larger  striation 
radii  decay  titm's  would  b(?  greater. 

f c*  \ 

We  have  investigated  drift-dissipative  modes'  ' (DD.M)  as 
a source  of  striation  diffusion.  The  DDM  whic;h  are  pictured  to 
grow  on  a background  provided  by  the  striation  structure  (.see 
Fig.  2)  .se;rve  to  diffu.se  the;  striation.  Our  results  ar<?  ba.sed 
on  the  assumption  that  DDM  are  not  sen.sitive  to  background 
electric  fields.  This  appears  likely  but  has  not  been  theoreti- 
cally verified.  The  fastest  growing  DDM  art?  typically  20  meters 
in  wavelength. 


To  understand  the  rang<>  of  applicability  of  DDM,  we 

define 


C 

s 


V 


k(T  + T,  ) 
e 1 

m, 


and 
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Figure  2 


St ri at  ion  from  Gradient- 
Drift  Instability 


DDM  Serving  to  Diffuse 
the  Striation 


G(?ometry  for  DDM  Diffusion. 


l-ll 


V . 
1 


V . + 0. 15 

in 


k2  c2 

y « 


Cl 


1 1 


where  k is  the  DDM  wave  vector  and  v. • is  the  ion-ion  collision 

y 


frequency.  Then  for  Li  »C  /v-  and 

1 ® 1 


» 2ti 


rsH  "'if 

\ ™e  V / ’ 


where  is  the  scale  length  for  density  variation  J_B , L||  is 
the  scale  length  for  density  variation  ||B,  and 

with  V the  electron-neut ral  collision  frequency  and  v^.  the 
electron-ion  collision  frequency.  Our  calculations  show  for 


k2  c2 
y s 


Cl 


= 1 


(which  we  are  free  to  choose)  and 


k = k \ 

z y \ w u)  , L,  j 
■^  \ ce  Cl  j ' 


for  the  DDM  growth  rate,  y, 


-1 

y ^ 0.07  sec  . 

‘-i 

8 6 

For  the  neutral  density  n -2.0  10  and  n^,  - 10  cm 

( 7) 

as  commonly  met  at  night  at  altitudes  above  -100  km'  the  con- 
ditions on  L|  and  L(|  hold  for  L|  « 1 km  and  Lj|  » 8 km.  Here 


-3 


1 

we  have  used. 


'll 


1 


V = 3.5  X 10~®  n + 10  ^ n^ 
e e 


-in  -7 

V.  = 8.3  X 10  n + 5 X 10  ' n 

X ^ 


T = T.  = 
e 1 


0 . 1 eV 


The  resulting  dir  fusion  coefficient  is^^^: 


D,  = 4?  = 0.07 


1 


1 ^ 


L|  = 0.07 


1 


Cl 


Here,  is  the  DDM  wave-vector  in  the  direction  of  density 


variation.  Our  calculations  show  k 


(p.  when;  is 


an  effective  ion-Larmor  radius  eciual  to  C /ij  ..  This  results 

s'  c 1 

in  a diffusion  time  estimate  of 

I.?  i-Tu.. 

_i  = .4  y ■ 

0.07^ 


For  a striation  width  ol  0.1  km,  corre.spondi  nu  tf) 
h?  U). 


Lj^  - 5 X 10  ^ km, 


c:  1 


0.07  C 


= 10  seconds 


This  corresponds  to  the  previou.sly  mentioned  e lec  t ron -i  on 
colli.sional  diffusion  tiriK'  of  about  1 hour. 

Put  another  way,  for  diffusion  due  to  electron-ion 
collisions 

,2 


V . c k ( T + T . ) 
D = ’ 

I 


v^.  C‘ 
ei  s 


ce 


q B 


(Jj  OJ  . 

ce  Cl 


where  is  the  electron-ion  collision  frequency.  Using 


ne  = 10  - Ti  = Te 


0.1  eV,  B = 0.5  gauss  results  in 


D 


'1 


8 X iqO  cm^/sec 


(9a) 


This  compare.s  with  the  drift-dissipative  estimate, 


D|  = 0.07  - 

1 


— = 2.5  X 10®  cm^/sec 
ci 


(9b) 
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under  the  sarm'  conditions  (with  n 1:  2.0  x 10  cm  '). 

3.2  rp;lation  between  the  size  of  various  parts  of  the 

STRIATION,  STRIATION  DECAY,  AND  DIFFUSION  COEFFICIENT 

To  hou r i St i ca 1 1 y treat  the  interplay  between  striation 

geoiTK'try  and  diffusion  coefficient  and  to  attempt  to  infer  the 

presence  of  DDM  we  consider  the  model  of  an  incompressible  flow 

► 

perpendicular  to  B to  which  a diffusive  contribution  will  be 
added.  Typically  striations  evolve  into  sheet-like  structures 
which  then  form  "bubbles"  or  break  up  into  parts.  Estimates 
apply  to  two  situations  which  appear  to  develop  in  striation 
evolution.  These  are  sheet  pinching-off  (Fig.  3a)  and  striation 
tip  steepening  (Fig.  3b). 

For  the  first  situation,  Fig.  3a,  we  consider  the  elon- 
gation of  a density  enhancement  in  the  y-direction,  due,  for 

instance,  to  an  ambient  electric  field,  E , in  the  negative 

ox 

,x-di  rect  ion  . We  take  the  magnetic  field  in  the  z-direction  and 

V • V = 0 (10) 

— ► 

where  v is  the  convective,  as  distinct  from  diffusive,  electron 
velocity  perpendicular  to  the  magnetic  field.  The  density  en- 
hancement has  characteristic  dimensions  L and  L at  t = 0. 

y 

Since  there  is  elongation  in  the  y-dir(?ction  and  since  the  flow 

-► 

perpendicular  to  B from  convection  is  incompressible, 

L^(t)  Ly(t)  = , (11) 

where  L (t)  and  L (t)  are  characteristic  dimensions  in  the  x 
X y 

and  y directions  at  time  t. 

Let  T be  the  time  when  the  convective  velocity  in  the 
x-direction  which  tends  to  make  the  striation  thinner  in  x is 
just  balanced  by  the  diffusive?  velocity  in  the  x-direction. 

Then 

3v  -9v  v., 

X _ y y 

9x  " 9y  Ly(T) 

■i 
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a)  Striation  Sheet  Pinching  Off 


b)  Striation  Tip  StoepeniiiR 


Figure  3.  Striation  Iso-Density  Contours 


where  is  the  velocity  of  the  striation  tip  with  respect  to 
the  striation  bast?  ( s«>e  Fig.  3a),  assumed  constant  in  time. 
Since  the  striation  distance  .scale  is  roughly  one-half  the 
striation  thickness 


li  ( T ) V ( T ) 

_x__  i:_  = V 

2 L (•  I ) L L y 
y X y 


But  from  diffusion, 


. - -1  r>  Sr*  n 2 

'■.X  n [ ) 


Combining  Eqs.  (12)  tuid  (13)  results  in 

,1/3 


VO  = 


The  corresponding  time  scale  is 

|2 


L 

X 


(0/2 


D 


1 


(4D|  ) 


1 /L  L 


( 12) 


( 13) 


(14) 


The  dependence  on  is  not  strong,  but  the  classical  time  scale 
from  Eq . (9a)  is  a factor  of  seven  greater  thfui  the  DDM  time 

scale  of  Eq . (9b),  ceteris  paribus. 

From  tnis  model  for  t i one  expects  L^(t)  = v^t. 

Density  decay  effects  can  be  estimated  very  crudely  by  noting 
that  provided  f ) Is  roughly  constant  (through  diffusion)  the 
ionization  enhancement  integrated  a magnetic  field  line  should 


decrease; 

as 

T/t 

t ^ T 

Using  D 1 

= 2.5 

10®  cm^/.sec, 

I,  = 4 X 10'’  cm,  L = 
X y 

-1  rv®  j 

10  cm  and 

< 

*< 

II 

o 

c;m/sec , 

values  which 

appear  appropriate  to 

the 
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( 9 ) 

oquat  or  i ill  sproad  F moasu  rorm'n  t s of  Kelley  et  al.  , we  obtain 
2 

r = 5 10  sec  tor  the  time  sc,ale  for  density  de.’cay  through 

q 

diffusion  and  L^( t ) /2  = 3 x 10  cm  for  the  corresponding  gradient 
scale  length.  A detailed  comfiarison  with  (?xpi>rimenl  al  data 
could  be  fruitful.  Mechanisms  such  as  striiition  splintering 
would  tend  to  decrease  r (see  Fig.  3c;);  the  effect  of  non- 
uniformity in  striation  thickncuss  as  a function  of  y ( sec'  Fig.  3a) 
is  ambiguous. 

The  second  situation  c.onsists  of  the  steepc'ning  of  the 
positive  y-tip  of  a striation  when  8E^/9y  0 where  F^  is  the 

x-directed  electric  field  at  the  striation  tip  (see  Fig.  3b). 

These  considerations  could  equally  well  apply  to  a "hole"  being 
driven  upward  by  gravity  in  equatorial  spread 

Provided  that  in  the  neighborhood  of  i>oints  A and  B 
(which  are  separated  by  a distance  Sy  along  the  y-axis)  is  gen- 
erated by  non-loc;al  .sourccis  it  is  reasontible  to  take' 

<5y  . (15) 

dt 


whore  n is  independent  of  6y  , since  the'  relative'  difference  in 
distance  f rcjm  a given  charge  source'  is  (Sy/b  where  L is  a typical 
non -steepened  cloud  dimension.  (From  the  one-dimt'nsi  on  al  ap- 
proximation as  studied  in  Volk  and  Ilaorende  1 ^ ^ ^ \ one  can  verify 
that  E^  is  continuous  acrcjss  a steep  dc'nsity  discontinuity  in 
y;  hence  that  it  is  not  generated  by  local  source's  in  y.) 


The  constant  of  priiport  i ona  1 i t y a is  givc'n  by 


a = 


3E 

c X 

B 3y 


(16) 


Now  let  y be  the  scale  lemgth  for  stt'epe'ii  i ng . The  diffusion 
ve^locity  is  Dj/6y.  Combining  this  with  Eejs.  (13)  and  (16)  give's: 

6y  B 9y 


6y 


(17) 


or 
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6d 


( 18) 


[i 


3E 

X 

9y 


Ti;  estimate'  c/B  (9E^/3y)  we-  approximate. 


c 

B 


B L 


where  (e/B)<5^  is  one-hal  f the  velocity  of  the  i nhomoftene  i ty 
with  respect  to  ambient  and  L is  the  characteristic  non- 
steepened  dimc'nsion  of  the  i nhomoftene  i ty . For  the  inhomoge- 
neity  rt'ferred  to  by  Costa  and  c<5  /B  = 8 x 10"^ 

cm/sec  and  L = 1 x 10  cm.  This  yields  (usinft  DD.M  diflu.sion) 


6y 


2.5 


10" 


10" 


10 


cm 


The  value'  from  electron-ion  colli.sional  diffusion  would  be  about 

d me.'ter.s.  Values  which  one  can  infer  from  Costa  and  Kelley's 

M 0 1 A A 

data^  ' are  in  the  ranR<‘  2 x lo  to  d x lo  cm. 
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DISCUSSION 


I 

1 . 

IX'nsity  i nhomof^ene  i t i e.s  whose  field  line  integrated  con- 
ductivities are  significantly  greater  than  ambient  (a  factor  of 
3,  at  least)  should  have  time  scales  for  decay  which  are  at 
least  a factor  of  sin"^  D (D  = magnetic  dip  angle)  lower  than 
those  for  the  homogeneous  nighttime  ionosphere  provided  the 
ambient  electric  field  and  neutral  wind  (as  a function  of  alti- 
tude; are  spatially  constant  and  substfin  t i al  ly  those  of  the 
ordinary  nighttime  ionosphere.  This  effect  is  caused  by  ioni- 
zation which  mtives  downward  and  is  not  due  to  di  f f us  ion  per- 
pendicular to  the  magnetic  field.  Even  in  the  absence  of 

downward  neutral  winds  the  characteristic  time  for  exponential 

2 

decay  would  hav'e  as  an  upper  bound  two  hours  times  sin  D. 

For  conductivities  not  significantly  greater  than  ambient  the 
time  of  two  hours  would  be  a more  conserv'at  i ve  estimate.  (We 
have  not  considered  heaved  ionization  which  could  return  to  the 
ionosphere  from  the  magnetosi)hero  on  time'  scales  of  several 
hou  r.s . ) 

; Decay  of  density  i nh<;mogene  i t ie.s  due  to  diffusion  per- 

■ 

I p«?ndicular  to  B is  strongly  entwine'd  with  convective^  non -di  I f us  i ve; 

( motion  and  in  addition  de-pends  on  whether  the;re  is  ce;l  1 isie)nal 

I diffusiejn  only  e;r  al.so  the  theore;!  i ca  1 1 y much  stronger  wave- 

particle;  scattering  diffusie;n.  The  latter  as  we  have*  noted  can 
be  due  to  dr  i f t-d  i.ss  i pa  t i ve  mode's  (DDM).  Altlunigh  pc;s.sibly  a 
! matter  e)f  obse;rvat  ional  interpretation,  observations  of  scale 

lengths  of  greater  than  100  meters  for  ionospheric  structur- 
ing^^*^’^^^  appe»ar  to  argue'  fe;r  a diffusie;n  me'chanism  which  is 
mucii  stre)nge*r  than  electron-ie)ii  collisional  diffusion. 
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Usage  of  a simi>le  modc'l  ( se«>  Fig.  3a)  which 
ation  stretching  linearly  in  time  in  one  dimension 
to  the  magnetic  fit'ld  yields  a characteristic  t im«' 
striation  decay  for  the  ionized  density  integrated 
magnetic  field  line  of 


entai Is  stri- 
pe rpendicular 
scale  for 
a 1 on  g a 


I 


L L 
X >• 


12/ 3 


(•in,  \1  V 


re  1 


(19) 


where  is  the  perpendicular  diffusion  coefficient,  is 

the  velocity  of  the  striation  tifi  with  respite!  to  the  striation 
base  and  and  L^.  are  t h«'  striation  dimt>nsions  at  t = 0.  For 
t » T,  providt'd  the  striation  dim<>nsion  perpendicular  to 
is  roughly  constant  as  a function  of  time  (through  diffusion) 
the  field-line  integrated  d*‘nsity  variation  i.s  at  r/t.  For 
n < 10^  cm  the  neuti-il  density  n <2.0  < 10^  cm  L =4.0 

^ , .1  o y 

km,  L = 10.0  km  and  v = 10  cm/sec,  t »=  3 x 10“  sec  with 

DDM  diffusion  and  t «=  1.0  » 10'^  sec  with  collisional  diffusion. 


Two  spatial  length  scales  appear  in  our  estimates.  The 
first  which  characterizes  the  x-dimension  perpendicular  to 
■s  t r i a t i on  e 1 on  ga  t i on  i -s 


Lx(t) 


^xS- 


re ' 


(14  ') 


with  variables  as  defined  in  Eq . (19).  The  second  which  char- 
acterizes the  steepening  at  the  tip  of  a striation  is 


6y  = 


( 18') 


Here  L is  the  characteristic  non-steepened  dimension  at  the  tip 
and  c£^/B  is  a nH-*asure  of  relative  velocity  within  the  striation 
prior  to  significant  steepening  which  is  taken  to  be  one-half 
the  velocity  of  the  striation  tip  with  respect  to  ambient. 
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Since  6y  corresponds  to  ;i  muximum  in  d<>nsi  Ly  <.>nhancement 
at  the  St  fiat  ion  tij)  one  would  exj)(?ct  variations  in  the  power 
density  spectrum  to  Ik;  associat('d  with  it.  If  the'ia;  are  mf)des 
at  rf'latively  short  wavelenpMhs  drivtin  by  lartte  density  f^radients 
(as  for  instance  DDM) , when  characteristically  6y  =0.1  km,  the 
value  of  E(j . (18')  c;ould  mark  an  uppc^r  bound  on  their  wave- 
lengths. On  the  other  liand,  if  Dj^  is  determinc*d  by  collisional 
effects,  when  characteristically  6y  = 1 meters,  the  value  of 
Eq . (18')  would  mark  a short  wavelength  cutoff  to  the  spectrum. 

It  is  i nt  ero;.st  ing  tf)  note  that  the  (vstimatt;  of  Eci . (18')  is 

very  similar  to  that  obtained  through  an  analysis  of  turbulent 

(13) 

behavif>r  by  Ott  and  Farley  ' although  the  physical  model.s 
appear  to  be  different. 

If  is  determined  by  collisional  effects,  the  structure 
of  the  striation  tij)  can  be  found  f>rovid(^d  w(?  assume  that  di  1- 
fusion  balances  convection.  Then 

Di  3n^^  _ - c y 

ii  9y  BL 

e 


D 


n^ , as  with  electron-ion  collisions  re.sults  in 


n ~ c 
e 1 


c <5  2 

X y_ 

BI,  2 


wh(»re  C.J  is  a constant.  If  Di  is  independent  of  n^^  , as  when 
collisions  with  neutrals  dominate. 


n^,  ~ exp  [-  c y^/(2BLI)|)J 

i.e.,  Gaus.sian  behavior  results.  have  not  found  the  wave 

spectrum  at  wavelengths  less  than  iSy  when  Di  is  determined  by 
turbu  lenc<> . 
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FUTURl-:  PRIORITIES 


1 


5. 

Work  in  the  I'ollowitu;  areas  would  be  useful; 

1.  Examination  of  nuclear  data  for  evidence  of  spread  F 
and  decay. 

2.  Calculation  of  decay  (through  falling  both  parallel 
and  perpendi cu 1 ar  to  B and  recombination)  for  a 
striation  whos»>  conductivity  is  large  compared  to 
the  background  for  disturbed  nucle'ar  c;onditions. 

n . Calculation  of  de<'ay  for  a .sti’iation  whose  conduc- 
tivity IS  not  large  compared  to  the  background  through 
falling  parallel  to  a field  line  under  disturbed 
nuclear  condi  ions. 

•1.  Veri  I K-at  ion  of  scale  lengths,  obse  rva  t i on  a 1 1 y , for 

e(iuatorial  .spread  F and  mi  d - 1 a t i t ude  spread  F.  Check- 
ing ol  pt“ rpeii d 1 cu  1 ar  di  I fusion  coefficients  through 
these  values. 

5.  Ik-terminat  i on  of  expe  r i m<‘n  l a 1 decay  behavior  for 

I equatorial  spread  F (which  seems  better  understood 

j than  mid-latitude  spread  F)  and  mid-latitude  .sprt'ad  F 

j and  determination  whether  stated  loss  mechanisms 

app 1 y . 

6.  Further  study  of  plasma  wave  phenormuia  — for  pf'rpen- 
dicular  diffusion  — and  determination  of  the  spectrum 

[ at  wavelengths  below  6y  when  Dj^  is  detc’rmined  by  wave 

ji  turbulence. 

1 


1-22 


'4 


I 

It  should  be  noted  that  barium  cloud  releases  at 
altitudes  of  200  km  would  tend  t<j  be  stable  to  DDM 
because  of  the  relatively  hlf^her  ion -neutral  c(jlli- 
sion  frequency.  Hence  decay  results  from  them  could 
be  misleading  if  applied  to  hif^her  altitudes. 


I 
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APPENDIX  A.  rONVI-:CTI  VI';  MOTION  PERPENDICPLAR  TO  B 
THROUGH  ELECTRIC  POLARIZATION 


U-siriK  a walerhapc  model 


(1)  . 


X = , p p 


(A-I) 


N = N, 


P P, 


N = , p > P, 


( A-2) 


where  p is  the  cylindrical  radius  around  the  axis  in  the  direc- 
* -► 

tion  of  B/|b| , from  Eq . (6)  we  have, 


Vj^  = 0 , P Pq 


V I i|;  = 0 , P > P, 


'1 


( A-3) 
(A-I) 


Using  e as  the  unit  radial  vector  and  e as  the  unit  azimuthal 
vector  , 


^ . Og 


Po-P 


= Vj^  ill  • e 


e 


c « p 


p +c 
o 


o 


(A-5) 


and 


E,  E„  - V,^  . e + N-  . 


' 1 

1 o 


1 U)  IB[  ' p 
Cl' 


P -L 

C) 


= I,,  E + 
2 o 


V X B 


— - 


+ N„ P rwi 

P ^ ‘^ci 


P</" 


( A-6) 


The  solutions  have  the  form, 
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z 


oo 

= ^ cos  ne  + p"  sin  iibJ  , ^ 

n=l  ^ 

Ejb  p ^ cos  n0  + b p sin  nG  , p > p (A-8) 
L n c ns  J o 


From  Eq . ( A-5 ) , 

^nc  ” \c  ’ *^ns  “ ^ns 

From  Eq . ( A-6 ) , 

/E  + V X B\ 

Ll  = (E,  - Z^)  [-^c )‘  ^ 

g X B • e 

+ (Nj-N2)e 

Define,  with  e^  the  unit  vector  in  direction  of  A, 

I »;  r , V X B ‘"l  ■ "2>  e J X B 
* ■ '‘'x  - * ^r-  * (Ij  - T^  IBT 


We  have, 


Z V,.p  • e - Z V|4-  • e = (Z^  - Z2)A  • e^ 

F ^ n -r  -L  D +e 


^Ic  ^ ^2  2 

i ^o 


cos  0 = (Z^  - Z2)A  cos  9 


Then 


(Z^  - Z2>A 

‘ic  = 


«^  = aic  X 
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A SYNOPSIS  OF  DRIFT  MODE  BEHAVIOR 
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1.  INTRODUCTION 


Drift  mod<-s  aro  ot  inttM-cst  in  HANPi  beh;ivi(jr  because 

they  furnisii  a mechanj  ,m  in  addition  to  the  usual  E 13  insta- 

bility^^^  which  can  allect  the  st  ructurintr  and  possible  decay 

of  hifihly  eiongatf'd  striations  (which  themselves  are  most 

likely  due  t(j  E ^ B modes).  Decay  of  E x B modes  into  drift 

—2  ( 2 ) 

modes  is  a iJossibD'  source  of  the  k power  density  spectrum. 

(k  i.s  modal  wave  number  in  direction  of  ob.se  rvat  i on  , which  we 

will  take  to  be  thf>  y-di  rec.'t  i on , B is  in  the  z-direction,  and 

( 3) 

Vn  is  in  the  x-d  i rec  t ion  . ) In  addition  it  has  bt^en  shown 
that  scale  lengths  of  0. 1-0.2  km  at  striation  tips  can  be 
attributed  to  drift  wave  turbulence.  (These  structures  are  be- 
lieved  responsible  for  anomalous  gigahertz  scintillation.  ) 

For  the  geometry  described  in  the  paragraph  above,  with 
E predominantly  in  the  x-di  rec.t  ion  , striations  tend  to  elongate 
in  the  y-direction.  E x B mode  grow't  h may  be  .strong  at  the 
striation  tips  where  one  expects  |1?  x Vn  | / | it  • Vn  j » 1.  How- 

J 

ever  along  the  elongated  striation  sides  one;  expects 
i E X Vn  I / I E • Vn  I «"  1 . It  i .s  then  h i gh  1 y likely  that  E x B in- 
stability is  suppressed  by  the  non-exi  stenix'  (jf  localized  E x B 
( 5 ) 

normal  modes.  (We  note  that  the  x and  y directions  in 

Kef.  5 are  intt'rchanged  with  thos'  used  here'.)  Drift  mtKlt^s 
provide  an  additional  instability  m«'chanism  which  i.s  not 
necessarily  .suppressed  wh('n  |E  x Vn|/|E  • Vn  | « 1.  Hence  it  is 
possible  that  the  resulting  plasma  turbulence  could  result  in 
the  decay  of  striations  eve*n  when  E x B turbulence  i.s  not 
indicated  . 

Within  thi.s  pa|)er  w<>  present  rt!sult.s  which  show  that 
drift  wave  instability  i.s  definitely  not  limited  by  E • Vn 
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A 


provided  Lc|E  • Vn|/(B|Vn;)  2 ' 10^  rm^/soc.  (Hero 

L = l(n  ^ dn/dx)  ^|.)  At  larpior  values  of  Lr  | E ‘ynj/CBlVnl) 
behavior  has  ru)t  yet  been  firmly  ost  ab  1 i sh<»d  . 

We  int('nd  for  this  synopsis  to  serve  as  an  introduction 
to  drift  waves  as  they  apply  to  striation  df^cay  in  the  iono- 
sphere and  as  used  in  Ref.  3.  We  first  present  a physical 
picture  of  neutrally  stable  drift  waves  and  then  indicate  the 
instability  mechanism.  We  next  show  the  connection  between 
t ho  CO  1 1 1 si  on  less  drift  instability  (universal  instability) 
and  the  dr  i f t-di  ssijiati  ve  instability.  For  a variety  of  iono- 
spheric conditions  we  then  pre.sent  drift-dissipative  modal 
growth  rates  and  frequencies  with  logarithmic  perpendicular 
der.sity  gradi en ts,  equa  1 to  (n  ^ dn/dx),  a factor  of  2.5  above 
the  minimum  logarithmic  perpendicular  density  gradients  produc- 
ing marginal  stability.  (.Mode  paramt'ters  and  k are  the  same 
as  at  marginal  stability;  modes  vary  as  exi)(-iwt  + + iky) 

Wf!  then  find  a simple  lOrm  for  a lower  bound  on  the  growth 
rate  for  drift-dissipative  instability  to  be  used  in  diffusion 
estimates.  Next  drift  dissipative  eigenmode  structure  is  in- 
vestigated for  use  in  diffusion  coefficients  and  the  origin  of 
the  diffusion  coefficient  D = 0.07  C“/w  used  in  Ref.  3 is 
indicated.  (C  = {k,(T.  + T )/m.'“,  k,  = Boltzmann’s  con.stant 
T = temperatun*  . m = mass,  i = ion  ciuant  ity,  e = elc*ctron 
quantity,  = cyclotron  frequency.)  Finally  the  possibility 
of  drift  wave  modification  when  E "Vn/IVnl  is  sufficiently 
large,  {(Lc|E  'Vnl/BlVnl)  S 2 < 10^  t^m^/secl.is  indicated,  and 
tin.*  relative  weakness  of  e li^ct  romagnet  i c effects  for  drift- 
dissipative  mode.s  a.s  contrasted  to  col  1 i si  on  1 ess  clrift  motles 
is  demonstrated. 
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liASIC  Ml'C’HANISM  FOH  NEUTRALLY  STAI^LK  DRIFT  WAVF:* 


I 


I 

r 


2.1  CHARACTHHISTICS  OF  BACKGROUND  PLASMA 

Fleet  ron  t emperatiiri'  is  , e](>ctrons  are  able  to  move 
alc-r.K  magnet  ie  field  lines,  backf^roiind  density  is  n,  density 
gradient  is  3n/3.x. 

2 . 2 MODK  STRUCTURK 

Modal  quant  itiees  are  indicated  by  subscript  "1",  vari- 
ati.jn  in  y and  z,  e U>ct  rostat  i c , (luasi-nout  ral  . density  (n^) 
and  pot.-nf  lal  1 e,x])(iky  + ik^z  - i wt ) . Flectrcjns  and  ions 

ari'  t r*  ••  d a^  f 1 ui  ds  . 

Therefore  v»>  1 oc  i t i 6’s  are: 

^’l  " ''"le  1 " ^’li  L 

or 


= - ik  c 1 


Lz 


k,  T 1 3n , 
h e 1_  ^ 


(with  Vp  the  sum  of  electron-neutral  and  electron-ion  collision 
frequencies,  ) . Diffusion  velocity  perpendicular 

to  l5  is  neplected  since  it  does  not  contribute  to  the  density 
equations,  Eqs.  (2)  and  ■.  V 

F’or  ions  with  no  motion  parallel  to  B, 


+ v.j^  • Vn  = 0 


We  note  j - 0 for  V E.j  =0  (e  l(?ct  rostat  i c modes). 


*The  di.scu.s.sion  of  thi.s  .section  and  of  Section  i.s  patterned 
aftc'r  Reference  6. 
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P’or  1‘U’c’trons 


3n 


1 


3t  +'’’-]'Vn  + nV*v^=0 
or  Kublrac’tinp  Kti  . (3)  from  Eq . (-1), 

V • V j az  = 0 

From  Eqs.  (2)  and  (5). 


(4) 


(5) 


q-t  J 

k.  t"  ~ IT 

h e 


(6) 


Substituting  E(| . (6)  into  Eq  . (3)  and  using  E<j . (1)  results  in 


— = _ — i<  T i — 
k ~ (jB  i)  (>  n 3x 


V 


De 


i.e.,  waves  move  at  eU'ctron  d i amagnt' t i c velocity,  , and 

are  neutrally  stable.  For  thermal  ec)  ui  1 i b r i uni  for  electrons, 

= pxp(-q^  ^l/kB^e)  . 

with  q^  (=  -q)  the  electron  charge.  For  q<l>^/kgT^^«  1 this  is 
the  same  as 

T rt^  ■ 

Hence  the  electrons  are  in  thermal  equilibrium. 

Further,  .since  n^  is  in  phase  with  and  E^  = -V4»j 
= -ikt).|.  Fig.  1 is  appropriate.  Referring  to  Fig.  1 the  flux 
in  the  x-di recti  on  is 

El  X B 

niC g — = ^ Inj  I Ui  j cos(ky  + k^z  - wt ) sin(ky  + k^z  - ust) 

Integrating  over  one  wavelength  in  the  y-direction  gives  zero, 
i.t?.,  no  net  transport  of  density. 
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3.  INSTABILITY  MECHANISM  FOR  DRIFT  WAVES 


Instability  of  drift  waves  f(jr  a fluid  plasma  can  he 
viewed  as  driven  by  deviations  in  the  ion  fluid  behavior  from 
that  described  in  Eqs . (1)  and  (3).  Drift  modes  are  unstaf)l6‘ 
when  >'  kv-  and  I w - 1 /kv,,  « 1.  (Here  m is  tin;  ro'al  part 

of  u)  (i*  + iy).)  Then  lags  n^  in  phase  along  the 

direction  of  wave  propagation  and  there  is  a mean  fluid  flux, 
corresponding  to  instability,  which  acts  to  decrease  the  back- 
ground density  gradient  by  transporting  fluid  from  higher  to 
lower  background  densities.  Hence  effects  which  act  to  reduce 
act  to  produce  instability.  Two  such  effects  are  ion 
inertia  and  finite  ion  Larmor  radius  corrections. 

On  defining  A/k  to  be  the  distance  that  4)^  lags  behind 
n.j  in  the  direction  of  wave  propagation,  the  condition  .sin  A > 0 
is  equivalent  to  a net  flux  against  the  background  density 
gradient.  We  first  show  that  y s kv„  sin  A,  then  that 

^ 0 corresponds  to  y ■ 0 (and  that  > 0 

correspond.s  t<j  y ''  0).  Finally  we  introduce  ion  behavior  and 
show  how  it  can  change'  w . 

3.1  EQUIVALENCE  OF  SIN  A > 0 AND  NET  FLU.X  AGAINST 

BACKGROUND  DENSITY  GRADIE.NT 

When  the  potential  lags  the  density  by  a distance  A/k 
in  the  direction  of  wave  propagation  the  geome>try  is  as  shown 
in  Fig.  2.  The  ptitential  is  represe'nted  as  4]^  = i'f’yl  <?xp(iky 
+ - iw^t  + iA)  and  the  density  is  n^  = exp(iky  + 

- iw^t).  On  taking  real  parts  of  K.j  = -V$^  and  n^  the  flux  is 

9 

cn.jEj  ''  B/B  = ]nj||4)jl  kc  {sin(ky  + k^z  - uij,t)  cos  A 

+ cos(ky  + k^z  - u)j,t)  sin  A}  c<xs(ky  + k^z  - w^t  ) . 
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On  avoraj>:i  nfj  over  wavelent^Ui,  w i t,h  ^ deno  t i n the  avera^'e  and 
I fixed,  this  results  in 

<(niVix)>=  |nj|4)^|  ke,  sin  A/2B  , (7) 

i.e.,  fluid  flux  in  the*  positive  x-direction,  against  the  back- 
ground density  gradient.  i Alternately  from  f'ig.  2,  remember- 
ing that  B is  in  the  positive  z direction,  ono  can  easily  see 
that  the  density  extremes  at  y = 0 and  y = tr/k  move  farther 
av-'ay  from  x = 0 when  A > 0.) 

3.2  RELATION  BETWEEN  y AND  SIN  A 

On  the  other  hand  if  the  ion  flow  is  given  to  first 
approximation  by  Ec} . (1)  and  f,  denotes  the  displacement  of  the 
ion  density  perturbation  from  x = 0 (see  Fig.  1 or  Fig.  2)  one 
has 

n^  = - f;  dn/dx 

and 

= clC/dt  . 

Y 1 

On  taking  C “ e cos(ky  + k^z  - u^t)  and  averaging  over  a 
period  in  time  one  has 

"l^lx  " ■ Y dn/dx)  . (8) 

On  using  Eq  . (0)  as  a first  approximation  for  j n 1 / | ({)  .j^  | and 
equating  Eqs.  (7)  and  (8)  this  results  in: 

Y = k Vj^^  s i n A . ( 9 ) 

3.3  EQUIVALENCE  OF  u - kv„  < 0 AND  |u)„  - kv,.,^|  « u),. 

WITH  instability'^ 

To  show  the  stated  equivalence  we  first  note  from 
Fqs.  (1).  (2)  and  (4) 
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-4) 


I 

f 

f 

I 

i 

i 

t 

I 

i 

P 


n 


1 

1 


(i  kv' 

Df 

i (A) 


( 10) 


2 

with  6 = k k.  T /ni  V . Usititt  Ix] . (6)  aRain  to  approximate 
z h <>'  e e 

Eq.  (9),  notinp  u)  " + i Y and  taking  the  imaginary 

part  of  Eq . (10)  results  in  (aftf'r  some  algebra), 


kvq,  iS(w  - kv^.  ) 
IX>  r De 


/u  ^2  2 

(kVj3^J  - “r 


(Y  + 


- kv„  6(w  - kv„  ) 

De  r De 

(Y  + 6)^ 


hence  the  equivalence. 


3.4  lO.N  BEHAVIOR 

Ion  behavior  within  drift  wave's  has  been  calculated  by 
Kennel  and  Hudson. Their  result  yields 

”l/'^  ^ ^De  - 

q(fj/kT^  w + RrpbCw  + i /b ) ’ ^ 

where  b = k^k^T^/ ( m.  , R,^  = T . /T^  and  = b(v.^  + 0.3bu..) 
= bv^  . ion-neutral  collision  trequency;  v^-  the 

ion-ion  collision  frequency.)  It  can  readily  be  verified  that 
Eq.s.  (10)  and  (11)  combine  to  give  Eq  . (lb)  of  Ref.  7.  On  the 
right-hand  side  of  Eq . (11),  the  contribution  -bw  in  the 
numerator  is  due  to  ion  inertia  and  the  contribution  linear  in 
R,pb  in  the  de>nominator  is  due  to  finite  ion  Larmor  radius 
effects.  For  large  electmn  conductivity  parallel  to  one 
ha.s  6 » 0) , kvjjj,.  Then  one  obtain.s  ( n ^/n ) /( c](j)^/kTp ) =:  1 and 
one  sees  that  both  ion  inertia  (-bw)  and  the  finite  ion  Larmor 
radius  term  (R,pba))  act  to  decrease  m , hence  to  cause  to  lag 
behind  n^  in  the  direction  of  wave  propagation,  hence  to  drive 
Instability.  Under  the  sanx'  conditions  it  is  easy  to  see  that 
collisional  contributions  in  are  stabilizing. 
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Wt‘  note  that  the  fluid  results  appear  somewhat  uncertain 
in  detail  in  that  in  calculations  which  we  carried  out  there 
were  di scrc^panc ies  from  Eq  . (11).  However,  our  discrepancies 
wt're  of  such  magnitude  and  sijtn  as  to  enhance  the  i n.s  tabi  1 i ty , 
so  that  the  results  presented  below,  which  are  ba.sed  on  Eq . (11), 
are  from  that  standpoint  at  worst  conservative  as  concerns  the 
effect  of  the  instability. 


I 


i 
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■1.  RKLATION  TO  OTHER  TYPES  OF  DRIFT  MODES 


Combination  of  Elqs.  (10)  and  (11)  is  also  in  agreement 
with  the  Vlasov  result  appropriate  to  each  species,  a, 
governed  by  the  equation, 


V X B • 


J_\f  = 
dvj  ^ la 


with  quas  i -neut  ral  i ty  , provided  with  v,p  = (kj^T^^/m  ) “ one  takes 
M ^ \j . » k V ')  » k \7  anH  u » /,i  The  inequality 


~ V.  » k^V,j,.  , 


V » k V'  , and  V » w . 
e z Te  e 


» 0)  corresponds  to  the  drift-dissipative  mode  for  which 
contributions  involve  the  plasma  dispersion  function  W in  the 
form 


) + i V 

e 

- no) 

ce 

/2  k 

z 

^Te 

index . 

The  C( 

V « (iJ  , 

V . « 1 

For  T = 0 . 1 eV  and  10'  < n,  <10  cm  we  have 

3 

\>  . = 1 .5  n /10‘  , and  tiie  limit  v » m seems  generally  reason- 
ei  e e 

able.  Tlie  limit  v » /2  k v~  is  the  isothermal  limit  of 

e z Tf' 

Hudson  and  Kennel  Comparison  of  their  Figs.  2 and  3 indi- 

cates agreement  iu-tween  t h<>  isothermal  limit  and  finite  heat 
conduction  drift  modes  within  a factor  of  about  3. 
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<1 


5. 


GROWTH  RATES  AND  PARAMETER  RANCJES  EOH 


DRIFT  dissipativt;  modes 


The  growth  rates  and  general  modal  behavior  stated 
below  are  subject  to  a number  of  uncf>rtai  nt  ies  involving  the 
precise  form  lor  electron  and  ion  behavior  as  well  as  usage 
of  fluid  results  for  b comparable  to  unity.  A major  Justifica- 
tion for  presenting  the  results  is  that  in  many  cases  the 

diffusion  coefficient  derived  from  them  can  be  much  larger 

( 3 ) 

than  that  resulting  from  coulomb  collisions. 

C'ombinaticjn  of  Ecjs  . (10)  and  (11)  and  solving  fcjr 
lu  ^ + iy  (d)^  = resonant  frequency,  y = growth  rate)  results 

i n 


and 


A + (A“^  + B 


U) 


r 


whf>  re 

T = 1 + x( 1 + y^) 


A 

B 


X 


2 2 2 
- Axy^  - (y/./2)‘^ 


2xyz  + (1  + y^)/2j 
, 2 


U)  . (i) 

c 1 ('e 


k2 


yz 


( 12) 


(13) 


y = kC  /w  . 

^ .s ' c 1 

/,  = C^/(Lvj) 
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Th*‘  condition  y = 0 relates  x,  y and  x.  at  marginal  stability. 
R(?KaftiinK  z as  a function  of  the  i nd('|)<,*nd('n  t variable's  x and  y 
and  minimizing  ■/.  succeuss i ve- 1 y with  respect  to  x and  y results 
in  three*  eeiuat  iein.s  for  x,  y and  z.  He>nce*  fe>r  give'n  ie:>nizf.'d 
de'nsity,  ne-utral  eie^nsity  and  i e)ne)sj)he>  r i c te'mpe'rat  ures  it  is 
pe)ssible*  te>  find  the-  lar^e-st  value  eif  L (=  L* ) at  which 
the're*  can  be  instability  and  the>  correspond  i nt?  values  ej  f x, 
y and  z. 

In  Table  1,  restrictinp:  wave'len^^ths  ale^ng  the*  maf^netic 
field  to  less  than  100  km,  fe)r  variable  ele'ctron  density  (NE) 
and  neutral  density  (NEUT)  with  L = 0.4L*  and  ~ ^ 

we  have  presented  y (RATE)  and  (FREQ)  for  values  of  x and 
y close  to  the>se?  which  resulted  in  L*  but  with  L = 0.4L*,  i .e.  , 
we  present  results  appropriate?  te?  marginally  stable  parame>ters 
but  with  spatial  gradients  a factor  e?  f 2.5  above  marginal 
stability  (LPERP  represents  L) . 

It  is  clear  that  there  are?  numerous  case.s  of  unstable 
behavieer  with  L on  the  scale  e? f tenths  of  a kilometer  v.'ith 
kL  » 1.  (We?  ne)te  that  the  conditiejn  kL  > 5 is  the  same  as 
I,y  > 2 < 10^  ( L in  cm)  . ) 

A simple  e>xpre'.ssion  fe?r  the  greiwth  rate?  obtains  for 
X = z z?  1 , with  y = 1.  One  finds 

y = 0.07  Cj^/L  . (14) 

The  result  is  significant  because  it  furnishes  a le?wer  bound 
on  growth  rates  which  can  be>  used  in  e.xpressions  for  the  diffu- 
sion coe*  f f ic  ien  t . The  value  y = 1 is  at  the  limit  of  validi'ty 
of  fluid  plasma  calculat  ie?ns  (y  < but  calculations  in  which 

the  ions  are  tre*ated  by  the'  Vlasov  equation,  col  1 i sie?r  lossly , 
but  with  y arbitrary  yield  a similar  re?sult. 
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6.  DIFFUSION  COEFFICIENT  FROM  DRI  FT-DI SS  I PATI V1-: 
MODES;  ADDITIONAL  PHYSICAL  EFFI'ICTS 


The  diffusion  coefficient  resulting  from  dri ft- 

2 

dissipaliv’e  mod('s  has  been  stated  to  be  rou<;hly  '(/^^  where 
is  a wave  numi)er  corresponding  to  the  modal  spatial  extent  in 
the  x-d  i rect  i on  . ^ ^ Hence  to  ('stimate  diffusion  it  was 

necessary  to  consider  the  eigenm(;de  equation  for  drift  modes. 
We  have  obtained  eigenmode  ecjuations  with  colli.sional  elec- 
trons ( x>  '^2’,'’Te’  €;ither  fluid  ions  or  co  1 1 i s i on  1 ess 

ions  with  no  restriction  to  y < 1.  Although  there  is  some 
uncertainty  as  to  the  precise  form  of  the  coefficients,  for 
both  cast's  thert'  is  an  eigenvalue  etjuation  of  the  general  form 


dx 


* -1 


with  8 of  ordf'r  unity,  (j  = ^ ^ ” ‘*'ci’ 

v^  = cE  • Vn/(  B I Vn  I ) . 

F’or  Lv  1 ''f  have  8 = Rfn(x),  L(x),  0|  and 

E(i  . (15)  ha.s  the  lorm 


-b  * i 
. 2 2 
dx  p j 


wi  th 


^ ^ 1 d^B  , ,2 

^ dx  ^ 2 ‘-2 

dx 


= 0 , 

( IB) 


^ 8 

dx  ~ L 


and 


,2  ,2 

dx  L 


Following  the  analysis  of  Rt*  f . 5,  st>e  t'specially  F.q  . (1),  we 
have  B'7p^  = a “ so  that  k ( a ( L(  . ) *.  Together  with 

1 X * 

Eq.  (14)  this  leads  to 
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n = >/k^  = 0.07  . 

,j  o 

^■sin^^  = 7 >c  10  cm/s('c  and  p.  = 3 x 10“"  cm,  one 

finds  ttiat  Lv^/(C^p.)  1 and  iience  tliat  (Electric  f luid  mf)difi- 

cation.s  art?  not  sif^nificant  for 

Lv  < 2.1  X 10^  cm^/sec 
c 

3 ‘1 

Thus  for  = 10'  cm/sec,  one  requir(?.s  L < 2.1  x 10  cm  for 

drift  modes  without  significant  plasma  velocity  shear  modi fica- 

1 3 

tion,  and  for  x ^ = 10  cm/sec,  one  requires  L <2.1  ^ 10'  cm 

for  the  same.  The  use  of  this  ineciuality  is  not  meant  to 
indicate  any  tendency  for  shear  to  iniiibit  drift  modes  when 
the  inequality  is  not  satisfied.  It  does  indicate  that  drift- 
mode turbulence  estimates  at  the  tips  of  striations  where 

« 0 art'  on  a stronger  theoretical  ground  than  those  along 
the  sidt's  of  striations  where  v / 0.  Locations  where  ~ 0 

are  thought  to  be  the  sources  of  anomalous  gigahertz  scintilla- 
( 4 ) 

tion.  Behav’ior  for  L\-  /(P;C^)  ^ 1 is  currently  under 

C IS 

i n vest i gat  ion . 

A final  point  concerns  the  limitatif>n  of  the  parametric- 

range  for  drift  mode  instability  due  to  e 1 c?ct  romagne  t i c coupling 

effc'cts.  Col  1 i sion  less  drift  modes  are  modi  fic'd  when  [1,  the 

ratio  ()i  plasma  pre.ssure  tc>  magnetic  fic'ld  prc'ssure,  is  greater 
( 9 1 

than  m^^/m.  . \V<?  find  for  dr  i f t -di  ssipat  iv(-  modes  that  there 

i.s  not  significant  c*  1 ec  t romagnet  i c modification  unless 
B > (m^v^J/(mjC^/L)  . 

E lect romagnc't  i c bc-havic^r  ari.se.s  bc-cause  ol'  the'  cont  ribu- 
tion of  the  modal  electron  current  along  ttie  background  mag- 
netic field  linos  to  the  modal  magnc't  ic  field.  The  electron 
current  varies  a.s  the  electron  fluid  velocity  parallel  to  . 
Hence  wc>  have, 


2-lH 


3v 


Ize 


(lE 


3t 


+ V V, 

o 1 ze 


Iz 


m 


Vo  i ^ 

m n 3 z 


or , 

qE,  k,  T . 3n . \ 

* Iz  b i‘  !_  1 J 

m ~ m n 3z  / ■ 

e e / 

Hence  the  mapnetic  field  source  term  is  a factor  w/v^  lower 
for  the  dri f t -diss ipat i ve  mode  than  it  would  be  in  the  absence 
of  electron  collisions,  and  it  plausible  that  the  threshold 
for  electromagnetic  modification  is  higher  by 
for  o)  = Cg/L. 


Ize 


( -iu  + V 


7 . SUMMARY 


In  this  synopsis  wo  havo  attempted  to  present  a basis 
for  the  physical  undorstandinp:  of  drift  modes  as  well  as  an 
indication  of  the  background  for  some  of  the  drift  mode  inputs 
to  Ref.  3.  For  C^/Lv.  » 1 and  kp^  ~ 1,  drift  modes  can  hav’e 
growth  rates  of  at  least  0.06  C^/L  = 0.1  comfiared  to 

cE^,/LB  for  E X B modes.  (We  recall  Vn  is  in  the  x-di  rect  i on . ) 
Further  dri  ft  modes  do  not  requi  re  E x }j  . Vn  >0  as  do  E x B 
nxjdes . If  one  neglects  electric  field  effects  in  the  direc- 
tion of  the  density  gradient  entirely,  from  Table  1 and  Ref.  3. 
drift  modes  .should  be  capable  of  preventing  .striations  from 
getting  thinner  than  0.1  km  in  many  case.s  of  interest.  The 
treatment  of  the  condition  E • Vn  ^ 0 (K^  0)  as  it  relates  to 

drift  modes  i.s  an  ongoing  problem. 
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1. 


INTRODUCTION 


I 

t 

i 

f 

) 

I 

I 


E X B modes,  analogous  to  the  fluid-dynamic  Rayleigh- 

Taylor  mode  have  long  been  recognized  as  a significant  element 

( 1-4 ) 

for  plasma  density  fluctuation  production  in  the  ionosphere. 
Observations  of  density  fluctuations  in  the  ionosphere  indicate 
some  tendency  for 


where  n..(k)  is  the  amplitude  of  the  one-dimensional  density 

^ (5-7) 

fluctuation  at  wav'e-number  k and  n is  the  background  density. 

Numerical  simulation  which  presumably  mimics  E x B mode  dynamics 

has  also  indicated  an  algebraic  variation  of  the  density  fluctu- 

ation  spectrum  with  k.^  ' A theoretical  explanation  for  the  k 


variation  based  on  the  decay  of  E x b modes  into  drift  modes 

( 9 ) 

has  recf'ntly  been  presented.  In  this  work  we  show  that  E x B 


tiKides  with  E /E 
x'  y 


non-zero  can  account  directly  for 


density  spectrum. 


We  consider  a background  (zero-order)  geometry  with 

spatial  variation  only  in  th.e  x-direction  (IB).  For  the  F region 

where  v.  /co  . 1 with  v.  th.e  ion-neutral  collision  freciuency 

in'  Cl  in 

and  0)  . the  ion  cyclotron  frequency,  the  plasma  velocity  is 

► 2 ^ I- 

V = cE  X B/B  while  the  current  density  is  j(x)  “ nE.  The  con- 
dition of  a stationary  background  results  in  3j  /8x  = 3/3.x(nE  ) 

= 0.  Hence  n = n(x)  and  non-zero  E^  at  any  x results  in  non- 
zero E = E (x)  for  all  x,  which  results  in  v = v (x)  or  a 

XX  y y 

velocity  shear  3v  /3x  = cE  /(BL).  [We  use  L = L(x) 

-1  -1  ^ ^ 

B (n  dn/dx)  ; i.e.,  L is  the  scale  length  associated  with 

the  amliient  density  gradient.] 
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I 
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I 

I 


( 


In  Section  2,  the  equation  for  the  E x B instability  in 
the  presence  of  plasma  velocity  shear,  obtained  by  Perkins  and 
Doles,  is  analyzed  generally  for  local  izf?d  mf)des.  in  Sec- 

tion 3,  with  k the  wave  number  of  the  E x B mode  in  the 
y-direction,  and  kL  »1,  the  equation  is  analyzed  for  the 
parameter  ranges 


(a) 


and 


E /E  « 1, 
X''  y 


kL  E /E  » 1 
x'  y 


(b) 


E /E  » 1, 
x'  y 


kL  E /E  » 
X'  y 


In  Section  4,  again  with  kL  » 1,  the  equation  is  analyzed  for 
the  parameter  range 

(c)  kL  E /E  ~ 1 

X y 


and  the  reason  for  the  erroneous  conclusion  of  Perkins  and 
Doles^^^^  that  E x B modes  are  stabilized  for  kD(E^/Ey)  > 2, 
kL  » 1,  is  shown  (D  is  a parameter  taken  roughly  equal  to  L)  . 

In  Section  5,  the  modal  results  are  summarized.  The  analysis 
for  range  (b)  for  nxjdes  of  the  form  of  Eq  . (4),  q.v.,  indicates 
if  the  ratio  E^/E^  becomes  too  large  (on  a scale  of  unity)  that 
spatially  localized  E x B modes  no  longer  exist.  Unstable 

modes  are  found  throughout  ranges  (a)  and  (c).  In  Section  6 

_2 

the  application  to  determination  of  the  k power  density  spec- 
trum is  presented.  The  eigenmode  structure  is  such  that  modes 
in  both  ranges  (a)  and  (c),  with  k large  on  the  scale  of 

E /(E  L),  have  the  same  dominant  variation, 

y X 


a k 


Mr>de  amplitudes  tend  to  be  limited  by  |dnj/dx|  = |Vn^|  because 
if  Idn^/dxl  > IVn^^l  modes  with  k*  » k will  tend  to  grow  more 
rapidly  than  (and  presumably  at  the  e.xpense  of)  the  mode  with 
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rij^(k) 

wi  th 
wh  ich 


_2 

This  leads  directly  to  a k spectruin 
the  process  of  decay  of  E x B modes  into 

_9  ( 9 ) 

can  also  lead  to  a k “ spectrum,  is 


Competition 
drift  riKJdes, 
readily  assesse 


d. 
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2. 


GENERAL  ANALYSIS 


Our  modal  analysis  is  based  on  the  model  of  Perkins  and 
Doles^^^^  which,  without  further  approximation,  yields  their 
Eq  . ( I'l ) which  has  the  form 


^ - ay  = 


(1) 


dx 


Here 


a = -k 


cE  B 


-1 


1 - 


n(  Y 


_ 1 yn"  ^ 1 / yn'  

Iw)  ~ 5 n(Y  - iw)  4yn(y  - 


(2) 


with  u)  = kcE^/B.  The  perturbations  in  plasma  density  n^  and 
electrostatic  potential  have  the  form 

, , , iky+Yt 

n^  = n^(x,k)  e 

, , ^ iky+yt 

= tl^^(x,k)  e 

X is  related  to  i|;j  by 

= [(y  - iw)/n]“  X 
and  the  prime  denotes  d/dx. 

We  first  solve  an  equation  of  the  form  of  Eq . (1)  under 
the  assumption  of  localized  modes.  Expanding  around  x = x^  one 
has 


a = + a'(x  - x^)  + ^ (x  - x^,)^  + . 


(3) 
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where  ) - a On  making  use  of  the  substitution  x - x = v 

U V J O 

^ seekin^t  modes  with  a sinpcl^-  peak  in  magnitude, 

Ixl  “ e"'^y  , Re  n '■  0 , one  readily  obtains  with  c a constant, 

•s  = -a  '/a  " and  : 


X = c exp  (a  ")^  i— 


. ^2  ^ 2a  ' 


.'2 

X ) + ■ „ 


J) 


(4) 


while  for  the  eigenvalue  condition: 


o 2a 


^ = 0 


(5) 


One  expects  x to  be  localized  within  a ranpe  x*  of  x 


such  that 


X - X 


2/2 


X* 


provided  that 


and 


Im(a" ) 


i, 


Im(a"  ) “ < 0 


Im 


7~^ 

(a  )- 


~ 0 


Condition  (6c)  can  be  written  more  precisely  as: 


llmia*")*!  x*x>  Im 


(a") 


J 


(6a) 


(6b) 


(6c) 


(6d) 


An  additional  check  on  validity  of  the  solution  procedure  is 
provided  by : 


3 

1 , , d‘  a . 

3 - \)>  T3 

dx 


« a 


(6e) 


since  higher  order  terms  are  neglected  in  K<]  . (3). 
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3. 


ASYMPTOTIC  ANALYSIS  FOR  kL  E ,/E  , » I,  kL  » 1 

y 


We  anticipate  that  modes  exist  for  which  y - iu)  ~ cEy/BL 

where  ~ denotes  equality  within  a numerical  factor  of  size  unity. 

Since  w = kcE  /B 
x' 


kL  » 1 
E 

y 


(7a) 


and 


Y = io) 


(7b) 


to  lowest  order  in  (kL  E /E  ) On  defining 

X y 

X = kx  ( 8 ) 

and  X = kx  , Eq.  (1)  becomes 

x"  ^-^X  ~ ® (9^) 


or : 


^10  ^ - ^o>  - ^o>^/2 


= 0 


(9b) 


where  = a^(x^),  prime  now  indicates  differentiation  with 

respect  to  x,  and 


, '"v  1 1 

y/w 

“ + 1 '■i"| 

1 p BLo)  (y/ui  - i)  4 

_kL(Y/u)  - i)_ 

2 ( y/u)  - i ) n 1 

(10) 


2 2 

Since  w/(y  - iw)  » 1 on  anticipating  n"/n  - 1/k  L the  final 
term  in  a-  can  be  neglected  compared  to  the  third  term;  further 
again  to  lowest  order  in  (kL  E /E  ) 

^ y 


Then  on  defining  a = Ey(  /E^(  , z =[kL(Y/w  - i)j~^  and 

using  Eq . ( 8 ) and 


so  that 


we  obtain 


du)  _ -u) 


(12a) 


dz  . 2 

— = - 1 z 

dx 


( 12b) 


- ( 1 - az 

+ z^/4) 

(13a) 

^ 2 
- 1 az  + 

i z^/2 

(13b) 

- 2 az^  + 

3 z'^/2 

( 13c) 

4 

6i  a z 

6 i z ^ 

( 13d) 

It  is  now  convenient  to  divide  the  analysis  into  the 
opposing  limits  of  a » 1 and  a 1.  We  note  that  it  is  simul- 
taneously possible  to  satisfy  kL  E^/E^  » 1 and  a»  1 provided 

kL  » a »1.  The  limit  kL  E /E  » 1 , (x  « 1 is  satisfied 

x'  y 

trivially  for  kL  *>  1. 


Limit  a » 1 


On  letting 


z z , 

a 2 
a 


(14) 


and  substituting  Ecjs.  (13)  and  (14)  into  (5)  (which  applies  to 
and  x as  well  as  a and  x)  we  obtain  on  collecting  terms 


Zo  = 4/3 


(15a) 


and  terms  ~ a 


Zj  = (4/3)^/^  a/|a 


(15b) 
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Since  is  real,  does  not  introduce  any  qualitatively 

different  ft'aturi's  into  the  eigenvalue  or  the  eiRenmfide.  The 
eigenvalue  is  given  to  lowest  order  in  a ^ by: 

^ = ^“o  ^ I Iff  - 

with  = w(x^^). 

On  using  Eqs.  (15a)  and  (15b)  in  (13b)-(13d)  one  has,  to 
lowest  order  in  a~^. 


a'  = - ia(z^ya)' 


(17a) 


a'"  = - 2a(z^/a)‘ 


(17b) 


= 6 i a ( z^ / a ) 


(17c) 


Hence  using  first  Eq . (4)  (with  in  place  of  a and  x in  place 
of  x)  and  then  Eq . (8),  y be  written. 


X = c,  exp  - ^ 7j 


1 ,4^3/2  k2(x  - 1^3  , _ 

2 3 M m ’ 


with  Cj  a constant . 

To  verify  that  the  rmjdes  are  consistently  localized  we 
note  from  Eqs.  (18)  or  (17b)  and  (6a)  that  one  has 

1/2  / 

x‘  = I ■ 


Hence  for  ( x - x 1 ~ x*,  contributions  linear  in  z^  (i.e.,  of 

^ — 1 ^ 

relative  size  |a|  compared  to  contributions  linear  in  z ) are 

of  order  a ^ compared  to  the  first  term  and  do  not  affect  the 

localization  provided  by  the  first  term.  Also  within  Eq . (6e) 


we  note : 


l(l/3)(x  - x^)  d^  a/dx^l 


2-  3 
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As  a further  check  we  note  that  the  localization  in  x imijlies 
localization  in  sinct; 

X*  9/9x  [(X  - ico)/n1"  ^ 

[(y  - iw)/n]'^ 

It  is  perhaps  relevant  to  turbulent  diffusion  that  the 
location  x^^  around  which  the  mode  is  "centered"  is  arbitrary 
and  not  determined  by  any  subsidiary  condition  (as  is  customary 
in  E X B wave  analysis)  such  as  d/dx(n~^  dn/dx)(x^^)  = 0.  Hence 
E < B i nstab i 1 i t i(;s  would  break  out  throughout  a re^^ion  of  non- 
zero density  gradient  of  appropriate  direction,  not  merely 
around  the  location  x,  althoufih  of  c;ourse  growth  would  be  most 
rapid  for  x^^  such  that  n ^ dn/dx  is  a maximum. 

Li  mi t a « 1 

Since  the  localization  of  the  modes  in  the  above  case 
depends  on  a x>  1 is  is  interesting  to  determini^  whether  the 
localization  persists  for  a « 1.  We  find  that  it  does  not  in 
general.  On  expanding 

z = + az^ 

in  Eqs.  (13),  substituting  in  Erj  . (5)  and  keeping  lowest  order 

9 

terms  one  obtains  z“  = -12/(6/s  + 1),  i.e.,  in  lowest  order  a'' 
is  positive  and  Im{(a..")‘'l  ~ |u|  at  the  large.st  . By  virtue  of 
Eq . (6a)  this  lt*ads  to  kx*  ~ |u|  “.  Then  on  using  Eqs.  (13)  in 

Ec] . (9b)  one  finds  that  the  successive  higher  order  derivative 
terms  in  (9b)  increase  at  least  as  rapidly  as  |a|~“  and  the 
localization  assumption  is  inconsistent. 

Hence  the'  calculations  in  the  two  limits  a » 1 and  a « 1 
indicate  that  localized  E « B instability  modes  are  associated 
with  sufficiently  large  values  of  a.  This  suggests  that  small 
values  of  u stabilize  the  Exp  instability  by  causing  the 


a 
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normal  modes  to  be  non-local  so  that  regions  of  unfavorable  as 
well  as  favorable  density  gradient  (see  Figure  1)  are  sampled 
by  the  same  mode. 


B 


Vn 

Favorable 
Density  Gradient 
for  E X B 
I n s t a b i 1 i t y 


Vn 

Un  f avorab  le 
Density  Gradi<‘nt 
for  E X B 
Instability 


dn 

dx 


= 0 


Figure  1.  Regions  of  Favorable  and  Unfavorable'  Density 
Gradie'nt  for  E x B Stability. 
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•1  . ASYMPTOTIC  ANALYSIS  FOR  kl.K  /L 

x'  y 


1 


kl,  • 1 


Fo  r this 
und  defiru?  B(  x 

' C) 


and 


with 


We  expect  S (x 
m o 

is  not  excluded 


case  we  anticipate  '/.^(x^^) 


.V,/(Y  - i 

) = Ey/(E^  kL)  . Then  within  E<i  . (9b)  we 


‘^l(Xo)  = - (1  - ^'1  ) 


r'l 


a{(Xo>  = - i 


a'(Xo) 


fx,,Zj  - 1 


s 


^2  = 

= 6/(kL) 

M 3/(kL)^ 

“5  = [^^2^^o^  ■ ■*] 

ug  = 26/(kL)^  . 

S (X  ) = (kL)"^  — ^ (-  — ) ■ 

">  dx"’-^  \"  dx/ 

1 in  Reneral;  the  possibility 


i . 

h a ve 
(20a) 

(20b) 

(20c) 

(21a) 

(21b) 

(21c) 

(21d) 

(21e) 

(21f) 

(22) 
= 0 
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The  expansion  for  is  clearly  valid  for  modes  with 
spread  Ax  «L.  Since 

= (kL)”^  (23a) 

iU)  d 

a'tx^,)  = (kL)'^  , (23b) 

from  Eqs.  (4),  (6c)  and  (6d)  one  expects  that  modes  should  be 
localized  around  the  region  specified  by 


Im 


(6c') 


with  eigenvalues  specified  by  . (3). 

From  Eqs.  (21a)  and  (2Ib)  we  have  to  order  (kL)  ^ in 
Eq.  (5), 


a 

o 


a /2a 


0 


(S') 


From  Eq . (6c') 


so  that  a > 0; 
o 


a'^/a  " > 0 
i . e . , Bz.j  > 1 or  , 


cE 

._Z  1 

BL  (y  - iWjj) 


> 1 


(6c '0 


For  cEy/(BL)  > 0,  this  imfilies  y - iw^^  < cEy/(BL).  However  it 

also  implies  y - iw  >0.  Thus  the  effect  of  the  velocity 

o 

shear  on  a localized  mode  in  this  region  is  to  lessen  the  growth 
rate  but  not  to  stabilize  the  instability,  and  the  question  of 
instability  (at  lea.st  for  localized  modes  in  the  form  of  Eq . (4)) 
becomt^s  (equivalent  to  th(*  e.xist(?nce  of  localized  modes. 


a 
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I 

k 

r 


With  a little  alRebra  w('  obtain  from  Kqs.  (20)  with 


H = (a'j^/2a 

",  (H 

' 0), 

H 

Zi 

i + 

22  2, 

U3Zi  - U2) 

2 

(«5Zii 

+ a. 

2 X 

6^'l  - ‘"4) 

For  012  = 0 

or  a_ 
5 

= 0, 

H cannot 

be 

rf'al  ; hence  Eq 

H 

= ^^1 

“2“3 

[^i  + (a^z 

2 

1 

a^)/(2a2U^z^)j 

"‘5 

[i  + (ug 

N 

1 

Posi t i vi ty 

( and 

hence 

reality) 

0 f 

H results  in; 

(24a) 


2 / “3 


3 “6 


^“2  “5 


CI2  “4 


2“3  “5 


(24b) 


(25a) 


But 


solving  for  requires 


z-  = 


1 “ “2''*3 


(25b) 


Eliminating  z^  from  Eqs . (25a)  and  (25b),  using  Eq . (21) 
and  writing 


6 


nl 

E Lk  nL 

X 


where  n,  and  L are  the  respective  values  of  n,  E^,  and  L at 
some  reference  point  in  x,  results  in 


IJy] 

1 

M 



^uL(  x^^)^ 

(3S2  - 1)^  (S2  + 1) 

^n(Xo)L^ 

1 (2S2  - 1)^  S2 

82(382  - 1)  - 2S3 

(26) 
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with  S2  evaluated  at  . The  value  of  S2  is  constrained  by  the 
condition  > 0 which,  combined  with  Eq . (25b),  results  in 

1 > a2a2/(a^a^)  > 0 . (26a) 

This  is  equivalent  to  S2(x^)  < 0 or  S2(x^^)  > 1/2. 

We  further  note 


cEy  (1  - 2S2) 

^ " ^“o  "bIT  (1  - 3S2) 


(27) 


It  is  clear  that  at  every  x^  with  the  right-hand  side 
of  Eq . (26)  positive  there  is  a single  "resonant"  k.  If  the 
right-hand  side  of  Eq . (26)  can  vary  continuously  for  positive 
values  from  0 to  infinity  subject  to  S,,(x^)  < 0 or  S2(x^^)  > 1/2 
then  localized  unstable  modes  are  in  principle  possible  at  all 
k.  Hence  the  existence?  of  regions  in  k without  localized  E x B 
modes  depends  on  the  form  of  n(x)  (which  determines  all  other 
variations  on  the  right-hand  side  of  Eq . (26)). 

We  now  survey  the  eigenmode  behavior  for  k first  large 
and  then  small  on  the  scale  of  (E  /EL). 

\ y/  X ' 


Large  k Limit 


Zeroes  of  the  right-hand  side  of  Eq . (26),  e.g. , 

^2^^o^  = -1,  locate  accumulation  points  for  modes  in  the  large 

k limit  within  the  ordering  kL  E /E  ~ 1.  We  consider  several 

X y 

functions  n(x)  to  determine  whether  such  large  k modes  exist. 

(For  such  modes,  with  S„  ~ -1,  y-iw  = (3/4)  cE  /(BL).) 

^ o y 

The  functions  are: 

n = a exp  (-  x^/2x^)  , Gaussian  (28a) 


n ^ 
n 


function  used 
and  Doles^l^l 
Eq.  (15)  with 
their  x 

o 


by  Perkins 
in  their 
X replacing 

(28b) 
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n = n exp(-Ax),  X = constant 


(28c) 


For  E(i . (28a),  S2  = -1  is  attainat)le  at  x = *x  with  x = -x 

corresponding  io  instability  for  > 0.  For  Eq . (28b),  S„  = -1 

^ 2 2 

is  attainable  for  x - x such  that  2d(x  - x)/D  = d - ( x - x)“/D  ), 
i.e.,  at  a value  (x  - x)/D  < 1.  For  Eq . (28c),  S2  = -1  is  not 
attainable  since  S,,  ' 0 . 

Our  results  for  this  limit  differ  frcjm  the  results  of 
Perkins  and  Doles^  which  indicate  stability  of  E x U modes 


for  E kD/E 
X ' y 


i.e.,  for  sufficiently  large  k within  the  order- 


ing E kL/E  -I.  (D  is  assumed  comparable  to  L.)  The  analysis 
^ ^ ( 10 ) 

of  Perkins  and  Doles  is  conducted  around  a point  for  which 

S,,  = 0;  at  such  a point,  appropriate  only  to  the  limit 

E kL/E  « 1,  there  is  nt)  mode  localisation  for  E kL/E  S 1. 

X y X y 2 

One  can  readily  see  this  f rc'm  the  inconsistency  of  (a')  /2a" 

real  and  positive,  Eq . (2-la),  and  = 0,  f 0 (from  S2(x^)  = 0). 

Further  we  note  on  defining  Eo  be  the  "approximate" 
solution  to  Eq . (20)  of  Perkins  and  Doles^^*^^  (and  hence  to  their 
Eq . (18),  that  if  one  assumes  |(X]^-X)/xl  « 1 that  x^  ~ X’ 
one  obtains  the  result 


d"(x-,  - X) 


,.2  cE 

D_  y 1 

'^2  BL  Tx  - 


whereas  for  self  consistency  one  would  expect  I « 1 
d are  as  defined  in  E(} . (28b).) 


(Here  D and 


Do  1 es 


(To  see  this,  we  note  that  in  the  notation  of  Perkins  and 
Xj^  satisfies  the  equation 


^ ^1  _ 1 kP(  Y*  - 1 + g"/4)  -t 


Y*  + a /: 
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I 


! 


I 

I 


On  usinp  their  Ec^s . (17)  and  (19)  with  = cE^/Bd  this 
becomes 


Yq  ig>^(x  - x^^)/d  ^ Y„(x  - 

Y - iw^  + (1/2)  wVay^cI^) 


(30) 


while  their  Eq . (18)  can  be  written; 


Yq  + ico^(x  - x^)/d  Y„(x 

Y - 


(31) 


We  now  replace  on  the  right  hand  side  of  Eq  . (30)  by  y ^nd 
check  the  consistency  of  x-j^  ~ X hy  subtracting  Eq . (31)  from 
Eq . (30),  thereby  obtaining 


d^(Xi 


X)/dx^ 


a2  2 

d x/dx 


On  taking  (x  - x^^) 


o^^^D^/(YQd^)  + ig)^(x  - x^)/d 

Y-ico  + iu)(x-x)/d 
' o o o' 


(D/k)^  and  D ~ d we  have  Eq . (29).) 


(32) 


Small  k Limit 

Poles  of  the  right  hand  sidf>  of  Ecj  . (26)  locate  rm^des 

with  k small  on  the  scale  E /(E  L),  i.e.,  kLE  /E  « 1 . The 

value  ~ d/dx(n”  dn/dx)  = 0,  reprt.'sents  the  standard 

no  shear  (E  = 0)  limit  with  Y - iw  = cE  /(BL)  and  a'/(a'0"  ►O. 

X _ o y / ' ' 

The  value  S„(x  ) = 1/2  corresponds  to  zero  growth  rate.  The 
condition  82(382  - 1)  - 28^  = 0 appears  to  present  the  possi- 
bility of  additional  localized  modes,  by  E(i . (27)  not  as 
rapidly  growing  as  those  at  82  = 0 . 
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5.  SUMMARY  OF  E x B MODES  WITH  BACKGROUND 
PLASMA  VT:L0CITY  SHEAR 

Effects  of  velocity  shear  on  E x B modes,  i .e . , of  a 

non-zero  value  for  E^,  the  backfjround  ele'ctric  field  in  the 

direction  of  the  i)ackground  density  Ki’^idient  , dejjend  on  the 

parameter  kLE^/Ey . In  all  cases  we  have  assumed  kL  » 1.  One 

may  distinguish  three  regime's:  kLE^/Ey  «1,  kLE^/Ey  ~ 1,  and 

kLE  /E  . » 1.  For  all  three  regimes  if  localized  modes  of  the 
X y 

form  of  Eq . (4)  exist  they  are  unstable  with  growth  rate 

such  that  0 < Y - cE  /BL. 

r y' 

The  regimej  kLE^/Ey  « 1 corresponds  to  the  standard  limit 

in  which  E is  neglected.  The  fastest  growing  modes  (othor.s 

^ -1 
may  possibly  exist)  are  localized  around  x with  d/dx(n  dn/dx) 

, -1  ” i 

= 0 with  Y = cEy/(BL)  + iou^  fuid  k^  = |n.j  dn^/dx|  ~ (k/L)“  « k. 

Modes  with  kLE^/Ey  - 1 are  determined  by  Eci  . (26)  subject 
to  the  constraints  that  the  right-hand  side  of  Ec) . (26)  is 

positive  and  either  ^2(x^)  < 0 or  ^^(x^)  > 1/2.  Corresponding 
to  each  value  x there  is  at  most  a single  unstable  mode  of  the 
form  of  E(} . (4).  Eigenvalues  are  given  by  Eq . (27).  For  this 
regime  it  is  noteworthy  that  k = [n7  dn.,/dx|  ~ k whenever  the 

A X i 1 

real  part  of  a' I {a")-  is  non -zero,  i .e  . , whenever 

Y - iu)^^  < cEy/BL.  Attainability  of  instability  at  a given  k 

is  a function  of  the  ambient  density  profile.  As  an  example, 

for  k becoming  large  on  the  scale  of  Ey/(E^L)  the  values  of  x^ 

approach  locations  specified  by  S2(x^^)  = -1.  Existence  of  this 

location  is  examined  in  Etps . (28a)-(28c)  for  three  commonly 

assumed  density  profiles.  The  limit  is  not  attain.*)^'?  for  the 

exponential  n = n exp(-Ax)  of  Eq . (28c)  but  it  exisi,s  for  a 

Gaussian  and  for  the  function  used  by  Perkins  and  Doles^^^^  in 

2 2 

their  Eq . (15).  It  is  attainalile  in  the  vicinity  of  d n/dx  = 0, 
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3 3 

provided  d' n/dx  ■ 0 which  is  generally  the  case  for  dn/dx  > 0; 

2 

and  d n / dx“  = 0 is  itself  attainable  whenever  the  density  has 

a maximum  and  is  ctjnstant  at  its  lowest  value  (at  "infinity," 

for  instance).  For  k becoming  .small  on  this  scale,  E /(E  L), 

y X 

the  behavior  for  the  limit  kLE^/E^  « 1 is  approached  directly 

at  S2  = 0 . In  addition  "small  k"  modes  appear  pos.sible  for 

.X  satisfying:  S„(x  ){3S„(x  ) - 1}  - 2S„(  x ) = 0.  The  con- 

''  ^ o ^ cy  oO 

elusion  of  Perkins  and  Doles^^®)  that  there  i.s  stabilization 
at  large  k ( in  this  regime)  is  due  to  an  invalid  approximation 
made  after  their  Eq . (20). 

The  limit  (a),  kLE  /E  » 1,  and  its  differentiation 

•k'  y 

from  the  limit  (b),  kLE  /E  ~ 1 but  getting  large,  can  best  be 

X y 

viewed  by  separate  consideration  of  the  cases  E^/E^  » 1 and 
E /E  - 1.  For  E /E  » 1 , the  dominant  terms  in  and 
derivatives  are  the  same  in  both  limits;  hence  the  dominant 
eigenmfJde -eigen value  description  as  given  in  Eqs . (16)  and  (18) 
is  the  same.  From  consideration  of  non-dominant  terms  in 

3 

Eqs.  (13)  and  (14),  limit  (a)  is  more  accurate  for  kL  > |Ey/E^|  . 

When  it  applies,  localized  modes  are  indicated  over  a wide  range 

of  x^^ , not  just  at  ^^(x^^)  = -1  (although  it  is  possible  that 

further  analysis  would  indicate  such  a tendency  for  limit  (b)). 

For  E /E  ^ 1,  limit  (b),  which  neglects  the  third  term  on  the 
y X 2 

right  hand  side  of  Eq  . (10)  and  hence  contributions  ~ 

is  inconsistent.  However  limit  (a)  with  E /E  « 1 which  retains 

> X 

the  asymptotically  dominent  terms  for  and  derivatives  is 
inconsistent  with  localized  modes  of  the  form  of  Eq.  (4). 

The.se  calculations  which  clearly  associate  values  of 
Ey/E^  large  compared  to  unity  with  localized,  unstable  E x B 
modes  sugg(?st  that  values  of  E^/E^  ^ 1 tend  to  suppress  insta- 
bility by  suppre.ssing  the  exi.stence  of  localized  eigenmodes. 
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(5.  RELI'VAI^CE  TO  TUP:  k~  POWKFi  DENSITY  SPECTRUM 


Unstable  behavior  is  not  seriously  inhibited  for  E /E 

> X 

large  compared  to  unity  and  one  can  generally  imagine  unstable 

E ^ B modes  with  very  large  valuers  of  kL,  i , kl.  » 1 whil<.* 

kp.  1 with  f the  ion  Larmor  radius.  {p.  = th  c. 

the  ion  thermal  speed  ( =/t . /m^  , with  T.  ion  temperature  and 

m.  ion  mas.s)  . } Hence  it  is  reasonable  to  assume  that  the 

parameter  range  kLE  /E  large  compared  to  unity  is  generally 

X y 

available  for  unstable  E x B modes.  For  this  range  within  the 
limits  of  either  Section  3 or  Section  1 the  dominant  eigenmode 
structure  is  the  same  and  is  given  by  Eq . (18),  in  particular 


|n^  dx  I ~ 3 ■ 

If  one  makes  the  assumption  (as  suggested  by  the  work  of 
Sleeper  and  Wei nstock ^ ^ ^ ^ when  applied  to  E x B modes)  that 
modal  growth  is  limited  whenever  the  modal  density  gradient  in 
the  x-direction  is  ctimparable  in  magnitude  to  the  background 
density  gradient 

I dn  j / dx  I ; n/L  (33 

_2 

one  obtain.s  a k power  dc'nsity  spectrum 


n^(k) 


k^L^ 


This  result  is  consistent  with  the  L dependence  for  a 
k pow€?r  spectrum,  ~ l/(k  L)  , indicated  by  Ott  and  Farley.  ' 
To  see  this,  on  defining  m as  the  running  integer  variable 
counting  the  number  of  m<ides,  we  note  that  the  contribution  per 
unit  range  in  k is  given  by 


n j ( k ) I dm/dk 
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(12) 

analot^ously  to  the  quantity  / fi(k,f3)  dO  of  Ott  and  Farley. 

If  one  introduces  no  nt^w  length  scales,  for  a (.slif^ht)  back- 
ground variation  in  the  y-direction,  one  has 

k « iTiTr/L  , m integral 


or 


dm/dk  a L 


from  which 

(n.j^(k)|^  dm/dk  = / f2(k,fl)  dO  « (k^L)  ^ 


The  criterion  that  modal  growth  i .s  limited  by  the?  condi- 
tion (32)  is  reasonable  from  the  standpoint  of  turbulent  cascade 
considerations,  since  ldn^(k)/d.xl  » n/L  allow.s  wav('  numbers 
k*  » k to  grow  more  rapidly  than  n^(k),  and  hence  pres\imably  to 

limit  the  growth  of  n^(k),  provided  the  value  of 

(E^  + E^^)/(Ey  + Ey;]^)  remains  small  whcui  |diiy/d.x|  ~ n/L.  {Here 

E , and  E , denote  the  sum  over  the  spectrum  in  k of  the  first 

X 1 y 1 

order  elec;tric  field  comi>on(M)ls  K^^(k,.x)  exp(iky)  and 

E ,(k,x)  <'.xp(ikv)  associated  with  n,(k).} 
y 1 • ‘■ 


In  thi.s  /nnection  one  note.s  that  the  mode  density  n^(k) 
and  e lec‘ t ro.stat  ic  potential,  , are  connected  by  Eij.  (9)  of 
Perkins  and  Doles,^^^^^ 


(y  - ik  cE^/13)ny  = ik  1);^  ( c/B)  dn/dx 

On  taking  y - ikcE„/B  = cE  ,/BL  and  k = [n"^  cltK/dx|,  when 

X y X I 1 

dn^/dx  « dn/dx  this  results  in 

ik  k^  U;y/Ey  = . (35) 


The  contributions  to  (K^j)"^  and  (Eyj)^  are  approximately  in  the 
ratio  1/3  for  kLE^/E^  ^ 1,  but  for  kLE^/Ey  «y  1 they  are  in  the 
ratio  k^/k(«l).  From  Eq . (35)  the  dominant  contributions  to 
and  (Ey^)^  occur  at  relatively  small  kL,  which  for 
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sufficiently  small  E ,/E  would  fall  in  the  ransf^  ELE  /E  « 1. 

X y X y 

Hence  provided  E,^/E^,  « 1,  the  ratio  (E^  + E^-j)/(Ey  + ) , 

which  could  turn  off  the  instability  if  it  became  sufficiently 
large,  tends  to  re?main  small  at  the  level  (j  f n^  at  which 
cascade  could  be  driven  unless  + Ey^  0. 

Eour  additional  points  may  be  of  interest; 

1.  Numerical  simulation  of  li  x B instabilities  in  a 

plasma  slab  as  summarized  in  the  five  cases  of  Table  I of 

Ref.  10  does  not  appe^ar  to  contradict  our  analysis.  Damping 

for  Case*  4 ( |E^/E^J  = 1.8)  is  consistent  with  the  breakdown 

of  mode  localization  for  relatively  large  |E,^/Ey|.  Faster 

growth  for  Case  5 (k  = 4tt/3,  |E  /E  I = 0.0)  than  for  Case  3 

X y 

(k  = 8tt/3,  |E,^/Ey!  = 0.6)  could  be  related  to  the  (relatively 

slight)  inhibitory  effect  of  shear  with  increasing  k.  (For 

kLE  /E  <K  1 one  has  cE  /BL  whereas  for  kLE  /E  » 1 one  has 
x'  y y'  x'  y 

3cE^./lBL.)  In  addition  the  greater  localization  and  hence 
difficulty  of  fully  accurate  representation  of  large  k modes 
with  a fixed  computational  mesh  spacing  could  cause  more  damp- 
ing at  larger  k. 

2.  If  one  ai)plies  the  c:rittirion  of  Ec}  . (33)  to  the 
limit  E =0  for  which  |n..^  dn.,/dx|  ~ (k/L)“,  one  obtains  a 

; i2-l 

, k spectrum  and  larger  values  for  |n^(k)/n|  . A k spectrum 

to  (jur  knowledge  is  not  observed.  On  th(>  other  hand,  the 

) limit  kLE  /E  » 1 would  .seem  in  general  attainable  for  E /E 

/ x'  y ^ x'  y 

small  compared  to  unity  and  L sufficiently  large. 

3.  The  estimate  of  Eq . (34)  can  be  compared  directly 

with  the  estimate  of  Chaturvedi  and  Kaw  based  on  limiting 
of  E X n mode  amplitudes  by  decay  into  drift-dissipative  modes. 
Detailed  calculations  indicate  that  the  threshold  for  drift- 
dissipative  instability  is  roughly  kn  .j  ( k )C . / ( n ) = 3,  with 
the  efff?ctiv«;  ion  collision  frecjuency  given  by  + 0 .3by^  , 

with  Vj^  the  ion-neutral  collision  frc-(iuency , ^ the  ion-ion 
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c-ullision  f rtMiucncv , and  b = k^C^/ui^. 

1 ' Cl 

this  leads' To 


n ^ ( k ) 


2 

JL  ^ 

2 2 

1 


Analogously  to  . (34) 


Comparison  of  Eqs . (34)  and  (36)  indicat(?s  that  F(i  . (34)  domi- 
natt's  whenever  L > (C^/v^)  ✓^/3. 

4.  At  relatively  long  wavelengths  for  which  kL  ^ 1 , 

since  our  calculations  are  all  for  the  limit  kL  » 1,  there  is 

no  particular  indication  that  the  power  density  spectrum  varies 
_2 

as  k . At  relatively  short  wavelengths  for  which  kp^  ~ 1 the 

fluid  picture  on  which  the  E x II  instability  is  based  becomes 

a poorer  approximation,  linear  stabilization  effects  are  more 

significant  and  the  availability  of  unstable  E x B modes  for 

fixed  L at  yet  larger  k to  limit  the  modal  density  gradient  is 

_2 

less  likely.  Hence,  here  too,  a k ^ spectrum  is  not  indicated 
by  our  arguments. 

_2 

To  summarize,  these  calculations  indicatt'  that  a k 

power  density  spectrum  for  k such  that  kL  » 1 » kp^  can  be  a 

natural  consequence  of  the  eigenmode  structure  for  E x B modes 

with  E^/Ey  f 0.  The  picturt^  is  that  E x b modes  over  a wide 

range  of  k are  unstable  off  the  pi-imary  gradient  length 

L = ( n~^  dn/dx)~  , that  they  grtjw  until  limited  by  a condi- 

_2 

tion  ldn^/dx|  “ n/L,  and  that  this  involves  a k power  den.sity 

spectrum  because  of  the  eigenmode  .structure  for  kLE  ,/K  , ~ 

X y 

i.e.,  Idrij^/dxl  “ kn^.  The  usage  of  the  criterion  |dn.j/dx|  = n/L 
is  of  course  approximate  and  more  detaiU'd  study  of  the  non- 
linear modal  f)ehavior  should  involve'  a fuller  accounting  of 
linear  properties,  the  dynamics  of  modal  interactions,  and 
other  possible*  turbulent  (or  wave  mode?  depe'ndent)  preice'sscs. 


'.I 
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